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SUMMARY 


Two methods of active pogo suppression (stabilization) for the space shuttle vehicle 
were studied analytically. The basis for both approaches was the linear quadratic reg- 
ulator, state space technique. The first approach minimized root-mean- square pump 
inlet pressure by using either full-state feedback, partial- state feedback, or output 
feedback with a Kalman filter. The second approach increased the modal damping 
associated with the critical structural modes by using either full-state feedback or 
reconstructed-state feedback. 

A number of implementable controls were found by both approaches. The designs 
were analyzed with respect to sensitivity, complexity, and controller energy require- 
ments, as well as controller performance. Practical controllers resulting from the 
two design approaches tended to use pressure and flow as feedback variables for the 
minimum-rms method and structural accelerations or velocities for the modal control 
method. Both approaches are suitable for the design of active pogo- suppression con- 
trollers . 


INTRODUCTION 

The problem of coupled feed-system - structural instabilities of booster vehicles 
(known as pogo) has persisted from early boosters to present vehicles. The classical 
approach to solving the problem has been to place accumulators on the suction side of 
the feed pumps. These accumulators, in general, have been tuned to be most absor- 
bent at the offending feed-system - structural frequency. Efforts have been made to 
expand the effective frequency range of such accumulators. As a generality, unless 
these devices are purely compliant (usually adding to existing pump inlet compliance), 
they introduce additional resonant modes to the feed system. These added modes can 
in turn couple with the structure to create new instabilities. 

In recent years, an attempt has been made to stabilize pogo by adding an active 
suppressor or controller to the feed system (ref. 1). By such an approach the sup- 
pressor response could be tailored to the required response over the bandwidth of the 
servo, and system flexibility could be improved. Farrell and Fenwick (ref. 1) began 
this approach by demonstrating an active suppression system experimentally on a test- 
stand, liquid-oxygen system. Then Lock and Rubin (ref. 2) considered the application 



of various suppression system designs to the space shuttle vehicle system. Seidel, 
Lorenzo, and Lehtinen (ref. 3) used the Lock- Rubin shuttle pogo model to study the 
problem of control- system design by using frequency-domain optimization techniques. 

In this report, time-domain optimization techniques and, again, the Lock-Rubin 
shuttle pogo model are used to design an active pogo suppressor. Two approaches 
are examined, both based on linear quadratic regulator theory: The first is called 
the minim um-rms design approach, and the second is called the modal control design 
approach. The objective in the minimum-rms design approach is to minimize the 
root-mean-square value of the high-pressure-oxidizer-pump (HPOP) suction pressure 
variations in response to a flow-rate disturbance at the HPOP inlet. Two variations 
of this approach were studied: (1) feedback of HPOP suction pressure through a 
Kalman filter and (2) partial-state feedback of selected variables (including HPOP suc- 
tion pressure) resulting in a suboptimal design. In the modal control design approach, 
the objective is to maximize the eigenvalue damping associated with critical structural 
pogo modes by feedback of reconstructed modal state information. The sensors re- 
quired for this purpose are a few selected vehicle accelerometers. 

The body of the report comprises three major sections. In the first, the develop- 
ment of a suitable system model and the controller design philosophy are discussed. 

A short review of linear quadratic regulator theory is presented next, foHowed by the 
second main section, which discusses the minimum-rms design approach. And finally, 
the modal control design approach is presented. An appendix describes the formulation 
of the space shuttle, state-variable pogo model, upon which the designs discussed in 
the report are based. 


POGO SYSTEM MODEL 

The space shuttle system configuration that was used for analysis and design of 
the active pogo suppressor is shown in figure 1. As in previous studies (ref. 4), only 
the liquid-oxygen side of the feed system is considered in the analysis. 


Overall Model Structure 

The liquid-oxygen tank is located forward in the main vehicle (fig. 1(a)) and is 
connected to the engines through a 31. 1-meter-long duct running alongside the main 
tank. This duct is represented in figure 1(b) by the portion between stations 1 and 2. 

(A more detailed schematic is shown in fig. 26.) The duct is followed by an elbow con- 
nected to an 8. 47 -meter-long crossover duct (stations 2 to 3) that feeds the low- 
pressure oxidizer pump (LPOP). The LPOP in turn feeds the HPOP through a short 
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interpump duct (between stations 5 and 7). It is at the inlet to the HPOP that the sup- 
pressor's piston actuator is assumed to be located. The HPOP exhausts the liquid 
oxygen through a short discharge duct into the engine thrust chamber. This simplified 
liquid-oxygen- system configuration is the same as the configuration used in references 
2 to 4. For instance, one major simplification is that the three shuttle main engines 
and their associated discharge ducts are modeled as a single engine of equivalent 
thrust . A more detailed development of this model is contained in appendix B . (Sym- 
bols are defined in appendix A. ) 

As in reference 3, the model used for this study contains seven liquid-oxygen- 
system structural elements that can move independently of one another. These ele- 
ments are the liquid-oxygen tank, the two feedline segments, the two elbows, the 
LPOP, and the interpump duct - HPOP - discharge line - thrust chamber combination. 
The propulsion system dynamics are characterized by flow rates and pressures 
throughout the system. 


Controls and Measurements 

Pogo oscillation is to be suppressed by means of an electrohydraulic piston actu- 
ator, at the HPOP inlet, that can insert or extract flow from the interpump duct so as 
to have an overall stabilizing effect on the system. This study assumes that measure- 
ments available to be fed back to the actuator are flows and pressures within the liquid- 
oxygen system and the accelerations of any of the aforementioned seven system struc- 
tural elements. 


State- Variable System Model 


Because a state-variable model of the pogo system had to be used in applying lin- 
ear quadratic regulator techniques, some modifications had to be made to the model 
discussed in reference 4. As shown in appendix B, the main modification was to for- 
mulate the line dynamics as a lumped- instead of a distributed-parameter model. 

Once this was done, the remainder of the task was simply to reformulate the set of 
ordinary differential equations into the standard state-variable form 


x = Ax + Bu + Dw 

(1) 

y = Hx 

(2) 


The details of this manipulation are given in appendix B . 
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Figure 2 shows schematically how the state-variable model is structured. Three 
main portions make up the overall model. Each portion has state variables associated 
with it. These state variables then comprise the overall system state vector x. 

The propulsion-sy stem-dynamics block in figure 2 contains flow and pressure 
relationships for the tank, pumps, lines, and engine. Coupled to it is the suppressor 
actuator, which is assumed to have first-order dynamics. Actuator input voltage u„ 

a 

will induce a change in actuator piston velocity and thereby produce a change in the in- 
jected flow rate Q & . The vector Xp shown in the figure is the propulsion system 
state vector and consists of the following 29 states: 

Xp = p^, q2 5 P2» Qgj P 3 » £ J4>P4» Qgs P 5> qg»P 6 > 


q 8 ,P 8 ,Q l ,P l’ ^2’ p 10’ q lV p ll> %’ P 4’ ^5’ P 7 ’ ^tc’ ^a) 

The capital Q's and P's are flows and pressures at stations indicated in figure 26. 
The lower-case q’s and p’s correspond to flows and pressures at discretely spaced 
points (lumps) along the first and second feedlines and reflect the lumped-parameter 
modeling process. This notation was chosen to be consistent with references 2 to 4. 
The variable Q a is the flow supplied by the actuator. 

The structural-dynamics block in figure 2 contains the modal structural dynamic 
relationships . A six-mode structural model was selected for this study . Although 
other studies (ref. 4) used one-mode-at-a-time analyses, we felt that it would be more 
advantageous to include the six predominant structural modes simultaneously in the 
model. In this way, control designs based on this model automatically take into ac- 
count any possible modal interactions. The state variables most convenient for analy- 
sis were chosen. They are the generalized modal displacements q and velocities q. 
Thus, the 12th-order vector of structural states is defined as 


T A 


= (w^'V 


v%) 


Because the vector x g (unlike Xp) is unmeasurable, it is important that the model in- 
clude structural variables that are measurable for use as feedback variables. The fol- 
lowing 12 measurable structural variables are defined as the structural output vector 

y s ; 



(*£!’ ^ 3 ’ 



Here, subscripts £1 and £2 refer to the first and second feedlines. For convenience, 
tank pressure P tb is considered as a structural output. 
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In judging the effectiveness of an active suppression system, and in designing con- 
trols that use reasonable amounts of control power, it is important to determine a suit- 
able model for external system disturbances. For this study, the main disturbance 
considered was a random noise signal that resulted in a flow-rate disturbance at 

the HPOP inlet . This disturbance is denoted by the signal w„ in figure 2 and is 
modeled by passing white noise through a 100-hertz, second-order filter with 0.707 
damping. (See appendix B for details. ) This modeling introduces two additional noise 
states, or * cn = <x onl , x cn2 ), where x cnl ^w cn . 

Thus, the overall 43rd-order state vector x is 

where the propulsion states x^ as well as the structural outputs y g are measurable . 
The elements in matrices A, B, D, and H are computed by subroutine POGMAT (ap- 
pendix C). 


Worst- Case-Model Parameter Selection 

Because the data (ref. 4) defining the model are preliminary and some model pa- 
rameters vary considerably over a flight, some method of handling parameter uncer- 
tainty had to be included in the control design process. We decided to select critical 
parameters so as to achieve a worst-case model. Worst case was defined as the model 
that was least stable in the open loop. The propulsion system parameters that most 
depend on flight condition are pump inlet bubble compliances and pump 

gains m^ + 1 and m 2 + 1. Ranges for these variables are shown in table I. Struc- 
tural system parameters that are required to define the model are mainly modal vector 
data, modal frequencies, masses, and damping ratios. These quantities are given in 
table II (also see ref. 4) for the two flight conditions of interest (least stable) - end 
burn and after solid-rocket separation. 

From preliminary stability (eigenvalue) calculations on the system A matrix, the 
sixth structural mode (27 . 2 Hz) was the least stable mode for end burn in the frequency 
range of primary interest (<30 Hz) for the overall coupled system. The 2.30-hertz 
mode was least stable for the after-separation condition. Therefore, the end-burn 
structural mode at 2 . 81 hertz (El in table II) was replaced by the 2. 30-hertz mode 
from the after- separation condition. The frequency of the sixth end-burn mode was 
lowered to 24 hertz, a condition that made for an unstable overall system mode at 
23. 8 hertz. Also, the pump gains were set to their maximum values and the bubble 
compliances were set to their minimum values (table I). 
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System Eigenvalues 

The eigenvalues of the A matrix for this worst-case model are shown in figure 3. 
The six modes primarily associated with the six structural modes of the uncoupled 
structural system are denoted by the solid symbols. The frequencies of these modes 
have not shifted appreciably from the original structural modal frequencies. However, 
the damping ratios for these modes have almost all decreased from their original value 
of 0.01. In fact, the mode at 23.8 hertz is unstable. Throughout this study, primary 
interest is focused on the 13 modes that occur at frequencies less than 30 hertz. There 
are seven somewhat lightly damped modes between 30 and 50 hertz that are associated 
with line dynamics . However, they are not affected by control or parameter variations 
and are included in the model mainly to provide better model accuracy at low frequen- 
cies. Three real poles are present: one for the actuator, one associated with the in- 
terpump duct, and one very high-frequency pole associated with the thrust chamber. 
Also present is the 100-hertz noise-coloring mode. 


System Frequency Responses 

To fully characterize the worst-case system model in the frequency domain would 

require a number of frequency-response plots equal to the number of system outputs, 

which could be as high as 41 . However, the stability characteristics can be observed 

by examining the plot of a key measurement variable, namely, HPOP inlet pressure 

. Figure 4 is a frequency-response plot of P ? to disturbance flow at the 

HPOP inlet. In this and all subsequent plots, pressures and flow rates are normalized 

and phase angles are usually principal values. The plotting resolution is 0. 1 hertz. 

Pressures are normalized to the nominal average value of thrust chamber pressure - 

taken as 2.07x10 pascals (3000 psi). Flows are normalized to the nominal average 

5 3 3 

value of total liquid-oxygen flow rate - taken as 3.52x10 cm /sec (21 500 in /sec). 
There are six main resonant peaks, all less than 25 hertz. The highest resonance oc- 
curs near the 2. 3-hertz, lightly damped, after-separation mode; and the second high- 
est occurs near the unstable, 23. 8-hertz mode. The presence of the unstable mode is 
evident in the phase angle plot, where there is a large phase lead near 24 hertz. 


Active-Pogo-Suppressor-Control Design Philosophy 

Figures 3 and 4 illustrate the pogo instability problem from two viewpoints - the 
modal viewpoint and the system- response viewpoint. In the modal viewpoint, demon- 
strated in figure 3, the primary means of evaluating stability is to examine modal 
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damping ratios. It has been found desirable in the past to have predominant system 
modal damping ratios all greater than 0. 01 in order to have acceptable overall system 
performance. This is not the case here. Five modes in the less-than-30-hertz region 
have damping ratios less than 0. 01, and of course one mode is unstable. The control 
design philosophy stemming from the modal viewpoint is this: Design a controller that 
will cause the damping ratios of all predominant modes to be greater than 0.01 with the 
restriction that the frequency of modes associated with actuator motion be less than 
some prescribed value. This restriction indirectly assures that the design be imple- 
mentable with an actuator that consumes a reasonable amount of power. 

In the system-response viewpoint, on the other hand, the primary concern is to 
minimize, in some sense, the response of critical system variables to outside distur- 
bances and at the same time to minimize the amount of actuator power required for 
control. Of course in minimizing system response, the control must, at least, stabil- 
ize the system. For the pogo problem, figure 4 illustrates a key system variable re- 

/v 

sponse to the main disturbance of interest, . The objective of a pogo suppressor 
designed from the system-response viewpoint is to reduce the resonant peaks of this 
response and to "quiet" the HPOP inlet pressure and, by doing so, to stabilize the 
overall system. This design approach requires that the actuator output response in 
the closed loop be computed to assess whether an excessive amount of control power 
is required. 

Both design viewpoints are discussed in the report. Although each has its advan- 
tages and disadvantages, either approach can be used to obtain a suitable suppressor 
design. 


LINEAR QUADRATIC REGULATOR THEORY 

Both approaches to suppressor design discussed in the previous section are imple- 
mented by using linear quadratic regulator (LQR) techniques. This section summarizes 
the main aspects of LQR and Kalman filter theory used in these design approaches. 

The linear, time-invariant plant to be controlled can be described by the state 
equations 

x = Ax + Bu + Dw (3) 

y = Hx + v (4) 


where w and v are white Gaussian plant and measurement noise vectors. The cor- 
relations for these noises are 

E{w(t)w(t + t) T } = ^ w <5(t) (5) 
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( 6 ) 


E{v(t)v(t + t) T } = ip y 6(T) 

where and ip v are power spectral density (PSD) matrices. 

The problem is to find a control u that m inimi zes a quadratic performance index 

J = e|x T R 1 x + u T R 2 u | (7) 

where R^ and R2 are, respectively, positive semidefinite and positive definite sym- 
metric weighting matrices. The solution to this problem is the so-called linear quad- 
ratic regulator (LQR) . If the state vector x is measurable, the resulting control law 
is 

u = -K x (8) 


where the control gain matrix K„ is 

v 

K c = < 9 > 

and S is the solution to the Riccati equation 

SA + A T S - SBR^ 1 B T S + R x = 0 (10) 

Figure 5 shows the structure for the L/QR solution with full- state feedback. The 
Riccati solution for obtaining feedback gains K„ can be obtained with existing com- 

V 

puter programs (ref. 5). 

If only the measurement y is available but not the complete state vector x, an 
optimal estimate x of the state is required. This estimate can be obtained by using a 
Kalman filter, described by 


x = Ax + Bu + K e (y - Hx) 


( 11 ) 


where the Kalman gain matrix is given as 

K e = PH T ^; 1 (12) 

and where P, the covariance of the estimation error e = x - x, is the solution to the 
Riccati equation 

AP + PA T - PH I V“ 1 HP + Di/> W D T = 0 (13) 
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Figure 6 shows the structure for the LQR solution when a Kalman filter is used in the 

feedback loop. Note that control gains K multiply x, instead of x as in full- state 

c 

feedback. 

Information on control system performance can be obtained by computing X, the 
state covariance matrix, a matrix whose diagonal elements are the mean-square values 
of the system states. When a Kalman filter is used, the covariance matrix X can be 
obtained by solving the Lyapunov equation 

(A - BK C )X + X(A - BK C ) T + (BK C )P + P(BK C ) T + D^ W D T - 0 (14) 

This equation can also be used with full- state feedback. For this case, P = 0 and the 

equation again can be solved for the covariance matrix X. 

Information on closed-loop system stability with full-state feedback can be obtained 

by getting the eigenvalues of A - BK . When the Kalman filter is used, stability is de- 

scribed by the eigenvalues of A - BK„ plus the eigenvalues of A-KH. These eigen- 

c c 

values are automatically calculated if a Riccati program of the type appearing in refer- 
ence 5 is used in computing the linear quadratic regulator and Kalman filter gains. 


MINIMUM ROOT-MEAN-SQUARE (rms) DESIGN 

The first of the LQR-based suppressor design methods to be discussed is desig- 
nated the "minimum- rms design" approach. As the name indicates, the stated objec- 
tive of this method is to minimize the rms value of the HPOP suction pressure . P ? to 
disturbance (flow-rate perturbation at the HPOP inlet). This objective can be 
achieved by minimizing mean-square values. Mathematically, the objective is repre- 
sented as 

Minimize = e|p|^ + r 2 u a) (15) 

Here, a penalty has been placed on commanded suppressor flow u a so as to limit the 
amount of control flow used by the suppressor. Standard performance-index weighting 
matrices R-^ and Rg can then be written as 

R x = diag(0, 0, . . . , 1; 0, . . . , 0) (16) 

26 

R 2 = r 2 ( 17 ) 
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Full-State Feedback 


So that baseline suppressor designs could be determined, a series of LQR solu- 
tions were run on the assumption that the full 43rd-order state vector was available 
for feedback. A number of runs were made with various values of r 2 . The computer 
program of reference 5 was used to solve the Riccati equation. For large values of 
r 2 , control gains K c and system eigenvalue damping ratios were small. For all 
values of r 2 the closed-loop system was stable, as predicted by theory. 

The open- and closed-loop behavior of one particular design case is shown in fig- 
ure 7(a). Here the open- and closed-loop frequency-response magnitudes for P ? 
responding to disturbance are compared. In addition to stabilizing the system, 

state feedback has greatly attenuated the frequency-response magnitude peaks in the 
frequency range 0 to 30 hertz. Nowhere does the closed- loop response exceed the 
open-loop response. 

The response of commanded suppressor flow u & to disturbance is shown in 

figure 7(b). Much of the control effort is, as would be expected, concentrated at the 

(•'o 

locations of the resonant peaks in the P^/Q^ response. The response amplitude is 
smallest at the points where minimums occur in the P ? response. Thus, in mini- 
mizing performance index J R , no real benefit is obtained by further reducing low- 
amplitude portions of the P„ response. The maximum in figure 7(b) occurs near the 
24-hertz peak in the P^/Q^ response, not near the 2.3-hertz peak. The reason is 
that the 24-hertz peak corresponds to an unstable open-loop mode but the 2. 3-hertz 
peak corresponds to a very lightly damped but stable mode. Thus, the control works 
hardest in mode stabilization. 

It is more informative to examine, for full- state feedback designs, the trade-off 

rsj 

between control effort and resultant rms values of P^ . For this evaluation, the mean- 
square system response was calculated for a number of designs. Each design corre- 
sponded to a different control weighting r 2 . All designs are optimal for any rms level 
of disturbance . In figure 8, results are presented for an assumed rms value for 


Q d7 ^Oq J of 10 percent. For larger allowable rms suppressor flows, the rms value 

of HPOP suction pressure decreases. Although not shown in the figure, ov, goes to 

*7 


infinity as rms suppressor flow is decreased to zero because the system is open-loop 
unstable. Hence without a stabilizing control, the rms value of P^ would be infinite. 

Another point to be noted in figure 8 is that the curve is asymptotic to a normal- 
ized suppressor flow of 10 percent, which is equal to the rms value of the disturbance. 
Thus, if a suppressor is available with a capacity approaching the rms disturbance 
flow, large attenuations in HPOP suction pressure can be achieved. This is, however, 
generaUy not the case. For example, a reasonably sized actuator may have a 15.2- 
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centimeter- (6-in. -) diameter piston with a maximum slewing rate of 254 cm/ sec 
(100 in/ sec), for a maximum flow rate of 46 900 cm 3 /sec (2860 in 3 /sec). If this is the 

O O 

3cr limit, the allowable rms suppressor flow would be 15 700 cm /sec (956 in /sec), 
or (normalized) 4.5 percent of the nominal liquid-oxygen flow. From figure 8 the re- 
sultant normalized rms suction pressure would be about 1.9 percent. 

Figure 8 represents an idealized situation, that is, full- state feedback and no 
measurement noise. Such designs are obviously impractical to implement without 
some simplification. One approach is to eliminate those feedback gains that least af- 
fect performance. This approach is not always possible; but, as in the section Partial- 
State-Feedback Minimum-rms Control Designs, reasonably simple state feedback de- 
signs can be achieved for the problem considered here. A second approach is to use 
only a few measurements of the most easily sensed variables and to estimate the un- 
measurable state variables by using a Kalman filter. This approach is discussed next. 


Kalman Filter Feedback 

Numerous measurement configurations can be used with a Kalman filter. Many 
considerations are involved in the choice of measurements - for example, sensor cost, 
accuracy, reliability, and noise levels. For the pogo problem, we decided to consider 
the simplest possible case, namely, a single sensor. The sensed variable chosen was 
P ? , a variable that has been shown in other investigations (refs. 3 and 4) to be a key 
feedback variable. 

A number of Kalman filters were designed that used P ? as a single noisy mea- 
surement. Each design corresponded to a different assumed level of measurement 
noise. The loop was closed around the filter by using the same sets of optimal control 
gains discussed in the previous section. The overall control loop is depicted in fig- 
ure 6. The results are plotted at discrete points and are connected by straight lines. 
The set of resultant design trade-offs can be displayed in an rms value plot, as shown 
in figure 9. This figure, as does figure 8, shows normalized rms HPOP suction pres- 
sure as a function of normalized rms suppressor flow for a normalized rms disturbance 
flow of 10 percent. The rms values were obtained by solving the Lyapunov equation (14). 
The curve from figure 8 is replotted in figure 9 and is now labeled "full- state feedback 
cases. " This curve serves as a lower bound to all the other cases, where a Kalman 
filter is used in the feedback loop. Three curves are shown, each for a different ratio 
of disturbance to measurement noise. These curves are also lines of constant Kalman 
filter gain. In all cases the normalized disturbance level is kept at 10 percent. The 

ru r>/ 

parameter for these curves is $ u ./$ tr , the ratio of normalized disturbance to measure- 
ment noise. As ip w / increases (less measurement noise), the curves approach the 


11 



full-state feedback curve. The dashed lines indicate loci of constant control gains; 
each locus is labeled with the control weighting r 2 used for the control design. For 
any fixed set of control gains, performance is sacrificed when a Kalman filter is intro- 

^ I /N/ rV 

duced into the loop: ov, increases along a constant-control-gain locus as ip Jib in- 

-t rj W V 

creases. 

One important difference between full- state and Kalman filter feedback can be 

noted by comparing the respective curves as ox gets large. For full-state feedback, 
~ -a 

ov, gets very small; but for Kalman filter feedback, the curves asymptotically ap- 

Y 

proach constant values of ov, . The reason is that the best a control with a Kalman 

7 

filter in the loop can do is to drive the estimate of P^ to zero. The irreducible error 

remaining is simply the error incurred in estimating . Shown in figure 9 are the 

asymptotic values of ov, for the three Kalman-filter-gain curves and one correspond- 

y 


ing to a ip^/ip^ of 1950 (curve not shown). These asymptotic values, which are ele- 
ments of the error covariance matrix P, are obtained automatically when solving the 
Kalman-filter Riccati equation. 

Figure 9 thus summarizes the rms performance of a large number of pogo con- 
troller designs; each employs P 7 measurement and a Kalman filter in the feedback 
loop. Kalman filter feedback is more complicated than full-state feedback because the 
measurement noise level must be determined. For example, suppose that the maxi- 
mum allowable value of Sq is 4 . 5 percent, the value used in the previous section on 
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full-state feedback design. Suppose also that sensor noise is such that ^ W /T v * is 19.5. 


And assume that the sensor noise is colored in the same way that the disturbance 

is colored. Then the rms value of sensor noise can be computed as 0.52 percent. 
Figure 9 shows that, for this design, normalized rms HPOP suction pressure would be 
about 3 percent, rather than the 1.7 percent for full-state feedback. For this case, 
the closed-loop signal -to-noise ratio on the P 7 sensor was 17 percent. 

Although acceptable performance can be achieved with Kalman filter feedback, the 
problem of complexity remains. As formulated, any of the Kalman designs requires a 
43rd-order system to be implemented. To determine whether simplifications are pos- 
sible, we are led to view the Kalman filter - feedback gain combination as a single- 
input, single-output compensator. 


Frequency-Domain Interpretation of Kalman Filter Feedback 

It is instructive to view the combination of a Kalman filter plus state-estimate 
feedback as a frequency-domain compensator (controller). The primary motive is to 
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see whether significant simplifications can be made in the controller. The following 
equations describe the controller: 

x = Ax + Bu + K 0 (y - Hx) (18) 

u = -K c x (19) 

where equation (18) is the Kalman filter equation (eq. (11)), repeated here for conve- 
nience, and equation (19) defines state-estimate feedback. Taking the Laplace trans- 
form of equation (18) for zero initial conditions and combining with equation (19) results 
in the transfer function matrix of the controller 

u(s) = -K c [si - (A - BK C - K e H)] _1 K e y(s) (20) 

The overall system, configured to show the feedback compensator, is shown in fig- 
ure 10. Here, u(s) = u a (s) and y(s) = P 7m (s) are scalars. Hence equation (20) can be 
expressed as 

V s) = G c (s)P 7m (s) < 21 > 

to define the controller transfer function G is). 

v 

The main variable parameters that define the characteristics of G (s) are the con- 
trol gain K and the Kalman filter gains K . These gains in turn are functions of r 2 
and ip w / ^ v , respectively. To examine the structure of G c (s), various controller 
transfer functions were calculated for various pairs of r 9 and ip /ip, . Figure 11 
compares three feedback-compensator transfer functions in terms of normalized 
frequency-response magnitude and phase. Commanded flow rate u a (s) and measured 
HPOP pressure P„ (s) were normalized as was done previously. The comparison is 

* 111 (V 

presented for controllers having the same set of Kalman filter gains (ip /ip = 1.95) but 
for three control weightings r 2 . All controllers are stable, as was determined by cal- 
culating the eigenvalues of A - BK - K H. This is a convenient result but is not al- 

L C 

ways the case with such Kalman filter - LQR-designed controllers. This fact was 
pointed out by Kwatny and Fink (ref. 6), who designed controllers based on a Kalman 
filter and found that the overall closed-loop system response was stable but that in 
some cases the controller by itself was unstable . 

The magnitude plot of figure 11(a) shows major peaks in the controller response 
near 2.3 and 23.8 hertz. This would be expected since this is where the least-stable 
system modes occur. The main trend shown is that, as control weighting r 2 was de- 
creased, controller gain increased - at all frequencies. Although not shown, similar 
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results were obtained when the 4 > w /4 ) v ratio was varied. In particular, when ^ w /^ y 
was increased, controller gain increased. 

The phase angle plot (fig. 11(b)) shows that the controllers all exhibit large 
amounts of phase lag. However, for smaller values of r 2 (higher gain controllers), 
phase lag is less. Controller poles and zeroes were calculated to determine why the 
controllers have such large amounts of phase lag. Each of the three controllers in 
figure 11 had right-half-plane zeros in its transfer function, which caused the addi- 
tional phase lag. Using a controller with right-half-plane zeroes is contrary to com- 
mon practice with classical frequency -domain techniques. However, Shaw (ref. 7) 
gives a root-locus example wherein right-half-plane zeroes are required to achieve a 
desired set of closed-loop poles. Thus, the existence of right-half-plane zeroes in the 
present results is unexpected but plausible. In particular, controllers designed for 

r 9 of 2x10 ^ and 2xl0 -4 both have 41 zeroes, 17 of which are in the right half-plane. 

" -5 

The controller for an r 2 of 2x10 has 41 zeroes, 13 in the right half-plane. The 
phase angles for these cases differ because the four additional right-haif-plane zeroes 
for r 2 of 2xl0 -3 and 2xl0 -4 occur near 24 and 40 hertz, causing additional 180° phase 
lags there (fig. 11(b)). 


Non-Minimum-Phase Controllers 

For a better understanding of the presence of right-half-plane zeroes, transfer 
function zeroes were calculated for various controller designs. Results for selected 
designs are given in table HI, which lists the number of right-half-plane controller 

r>j rv/ 

zeroes for various combinations of r 9 and ib / ip . The conclusion that is suggested 
by the table is that, in general, the higher the controller gain, the fewer the number of 
right-half-plane zeroes. Higher gain means either higher optimal control gains or 
higher filter gains . Further down the table and to the right, the number of right-half- 
plane zeroes decreases. In fact, the highest gain controller = 195, r 2 = 4x10"®) 

has no right-half-plane zeroes at all. From these results, we concluded that a low- 
control- effort stabilizing pogo controller will tend to have right-half-plane zeroes. 

This somewhat counterintuitive conclusion was also reached by WidnaU (ref. 8) in a 
sampled-data optimal control design of the Apollo autopilot. There are similarities 
between his problem and the pogo stabilization problem in that both systems to be con- 
trolled are unstable and contain a number of lightly damped modes. Thus low-control^ 
effort optimal control of unstable plants appears to lead to non-minimum-phase con- 
trollers. 

The remaining problem was to simplify the controller designs to a transfer func- 
tion of reasonable order (10th to 15th order). A method developed by Hutton and 
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Friedland (ref. 9) was used on a representative transfer function for an r 9 of 2x10 
and a ^ w /^ v of 1.95. The technique produces a series of reduced-order transfer 
function approximations, beginning with 43rd order and extending down to any arbitrary 
lower order approximation. For reduced transfer function orders of less than about 30 
the transfer function fits were poor, probably because the non-minimum-phase charac- 
ter of the transfer function made simplifying difficult. Thus, we decided not to pursue 
the Kalman filter approach further but instead to look into simplifying the full- state 
feedback designs by eliminating certain of the feedback variables. 


Partial-State-Feedback Minimum-rms Control Designs 

Sophisticated output-feedback techniques for simplifying LQR designs were dis- 
carded in favor of a much simpler, intuitive approach. Gain simplification was per- 
formed as follows: The feedback variables consist of four types - flow rates, pres- 
sures, structural velocities, and structural displacements. Thus, gains were grouped 
into four classes, and plots were made of each class for various LQR designs. In all 
plots, gain values were normalized. Velocity gains were normalized by the acoustic 
velocity in liquid oxygen (51 800 cm/sec; 20 400 in/sec), and displacement gains were 
normalized by the line lump length (345 cm; 136 in.). Gain (absolute) magnitudes were 
compared within each group. Then only the gains with the largest magnitudes were re- 
tained for use in the simplified gain control. It is difficult to rationally compare gains 
from one class with those from another. As discussed later in this section, other 
rationale were used to eliminate complete classes of gains. 

Figure 12 summarizes the results. Figure 12(a) shows the flow-rate gains as a 
function of the inverse of the control weighting parameter r 2 - Here, gains for states 
25 (LPOP exit flow rate), 27 (thrust chamber flow rate), and 29 (suppressor output 
flow rate) are obviously predominant. Thus, we decided to discard the rest of the 
flow-rate gains and keep these three. Figure 12(b) shows normalized pressure gains 
as a function of r 2 ^. Most gains increase with decreasing r 2 , but the most notable 
result is that gains for states 24 (LPOP suction pressure) and 26 (HPOP suction pres- 
sure) are dominant. This result is not unexpected because HPOP suction pressure was 
weighted in the performance index. Figure 12(c) shows the gains for the generalized 
displacements. The dominant gains correspond to the displacements of the last two 
modes. From figure 12(d), it is not apparent that any of the normalized, generalized 
velocity gains are dominant. For this reason and so that measuring (or reconstruct- 
ing) all the structural velocity states could be avoided, we decided to eliminate all 
structural variables from the feedback loops. There remain two gains unaccounted 
for: those associated with the noise-coloring states. These states are unmeasurable 
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and to include them would have meant implementing an observer to reconstruct them; 
thus they too were dropped. 

The gains left then were those for LPOP inlet flow and pressure, HPOP inlet 

pressure and outlet flow, and suppressor outlet flow. The performances of these 

simplified-gain, minimum- rms control designs are compared with those of the full- 

state feedback designs in figure 13. The rms performances of four simplified gain 

designs are compared with the performances of the full-state feedback designs upon 

-3 -4 -5 

which they were based. The base r 0 values chosen were 2 ><L0 , 2X10 , 2XL0 , 

—6 " 

and 2x10 . All simplified gain designs were stable, a fact not guaranteed by the gain 

elimination process. For any fixed value of ox , the simplified gain control is infe- 

v^a 

rior to the corresponding full- state feedback control. Controls compare least favor- 
ably at the higher values of rms suppressor flow. However, for the lower values, for 
instance for a ox of 4, the simplified design is not much worse than the full-state 

feedback design. It is in this region of that the practical upper lim it on o^ 

would probably lie. Thus, it appears that simplified gain designs are a reasonable 
solution to the suppressor design problem as long as the lower values of o^ are 
used. a 

As a check on the results shown in figure 13, a closed-loop frequency response 
was run for the lowest-gain simplified gain control shown in figure 13 (based on the 

_3 

full-state feedback control designed for an r 2 of 2X10 ), whose normalized rms sup- 

pressor flow requirement is about 4. Figure 14 compares this response to the open- 
loop frequency response. Closed-loop performance is quite acceptable because the 
disturbance is everywhere attenuated to values below open loop, as was the case for 
full-state feedback. Comparing figures 14 and 7(a) shows that, except for the region 
near 25 hertz, the simplified feedback design is very close in performance to the full- 
state feedback design, and both designs use essentially the same amount of suppressor 
flow. 

A major concern with any control system, optimal or otherwise, is how sensitive 
the closed-loop response is to changes in system parameters that may be imprecisely 
known. Therefore, we decided to check the sensitivity of the closed- loop response of 
the simplified gain design shown in figure 14 to changes in key system parameters. 
Experience has shown that two parameters that have major effects on system stability 
are the frequencies of the first and sixth structural modes. Thus, the frequencies of 
these modes were varied, and the frequency responses were calculated for the system 
that uses the simplified gain controller. Figure 15 shows a check case where the 
frequency of the sixth structural mode (23.8 Hz) was increased by 10 percent. There 
was no serious performance degradation. Similar results were observed when the 
first structural modal frequency was varied. 
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In summary then, it has been demonstrated that a partial- state feedback control 
for pogo suppression can be designed. The gain elimination procedure used was effec- 
tive in eliminating less important feedback variables. For higher performance con- 
trols, those that require large amounts of control effort, the simplified controls suf- 
fered appreciable performance degradation. However, for lower performance con- 
trols, simplified gain control was nearly as good as full-state feedback. Sensitivity 
of the simplified controls does not appear to be a significant problem. 


MODAL CONTROL APPROACH TO SUPPRESSOR CONTROL DESIGN 

The second approach taken in designing a control to suppress pogo oscillations is 
designated the "modal control design" approach. As already mentioned, the aim in the 
modal control approach is to increase the damping on the least stable or most sensi- 
tive system modes while not unreasonably increasing the modal frequencies associated 
with the actuator. The focus in this design method is primarily on the degree of sta- 
bility, with secondary focus on minimizing system response to disturbances. 

The LQR and Kalman filter techniques lend themselves directly to this approach. 
The key difference between the modal control approach and the minimum-rms approach 
is in the selection of the quadratic performance index to be minimized. The difference 
in the character of the index then leads to differences in the resulting control systems. 


Performance Index Formulation 

The first step in the modal control approach is to obtain the modal form of the 
system state-variable model. For the system state-variable equation 

x — Ax + Bu (22) 


a modal state vector 17 can be defined as 

rj = T -1 x (23) 

where the matrix T - ^ effects the transformation from the original system coordinates 
into modal coordinates. The modal state equations can be obtained by substituting for 
x in equation (22) from equation (23): • 

r? = An + T -1 Bu (24) 


where A = T -1 AT. 
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One convenient form of modal state equation can be obtained if the matrix A is 
the block-diagonal form of A, whose structure is 


An 


0 


A = 


( 25 ) 


0 


A,, 


An element A^ is either (1) a scalar equal to a real eigenvalue on of A or (2) a two- 
by-two matrix of the form 



(26) 


where a . and /3. are real and imaginary parts of a complex eigenvalue pair of A. 
With A being block diagonal, the matrix T becomes the modified eigenvector matrix 
(ref. 5). The block-diagonal, modified eigenvector matrix form of the eigenvalue- 
eigenvector problem is used here for computational convenience because no complex 
arithmetic is required. 

The quadratic performance index that allows the designer to obtain a control that 
tends to increase the damping on certain selected system modes is set up as follows: 
First, the eigenvalues of the matrix A are obtained. From the eigenvalues, a block- 
diagonal matrix A is formed. Next, the mode whose damping is to be increased is 
selected. There are then two modal states and V£ + i> corresponding to this mode, 
that are to be altered. Penalizing these states in the performance index will cause 
their amplitudes to be attenuated. As it turns out (especially for lightly damped modes) 
an attenuation in and t^ + 1 will lead to an increase in damping of the associated 
mode. Thus, the "modal" performance index can be written as 


j m = fj" [ 1 i ( v e + ’’j+i) + uTR 2 u ] dt 


(27) 


where the integral form is used for since random disturbances are not considered 
explicitly in the modal approach. Now the remaining step is to express and 7^ +1 
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in terms of the original state vector x. It is assumed that matrix T _1 has been cal- 
culated. Define 



(28) 


A 



where tj are row vectors. Then = t^x and rj £+1 = + ^x, from equation (23). 

Now if Tfy and Vg + i are substituted into equation (27), the standard form 


J M = f” ( xTR l x + uTB 2 u ) dt 
is obtained for J M> where the elements of R.^ are 


(29) 


( r l ) i>3 " + Wl.iVl.j) (30) 

To control the damping on other modes in addition to these, simply identify the indices 
of the associated modal states (m, m + 1, p, p + 1, etc. ) so that the elements of R^ 
become 



(31) 


where k takes on the values of the first indices of penalized modal state pairs. 
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Full-State Feedback 


As already discussed, the two least stable modes of primary interest in the liquid- 
oxygen-system model were those at 2. 3 and 23.8 hertz. These modes are referred to 
as "critical modes" (fig. 3). Other modes, which although lightly damped were not 
considered, were those at 4.72, 22.5, and 8.38 hertz. These modes are very insensi- 
tive to control and system parameter variations and never become unstable . Thus, all 
effort was concentrated on increasing damping on the 2.3- and 23.8-hertz modes. 
Full-state feedback designs were calculated for the case where both critical modes 
were equally weighted. The primary criterion for judging effectiveness is the increase 
in modal damping on these two modes. 

Table IV shows the results obtained for two different control weightings r 2 - For 
a weighting of 0.1, damping on both modes was in excess of 10 percent, which is gen- 
erally considered an acceptable target design value. All other modes were essentially 
unchanged. For both cases, the 50-hertz actuator pole was not moved appreciably. 
This shows that the two modal designs quite efficiently concentrate available control 
energy on only the two selected modes (2. 3 and 23.8 Hz). 

Eigenvectors associated with these designs were plotted so that insight could be 
gained into the behavior of the modal designs. The propulsion system portions of the 
eigenvectors (corresponding to the states that are either pressures or flow rates) were 
plotted as either pressure mode (figs. 16(a) and (b)) or flow mode (figs. 16(c) and (d)) 
shapes. The mode shapes are simply the normalized magnitudes of the (pressure or 
flow rate) eigenvector components plotted as a function of the spatial location of the 
corresponding pressure or flow-rate variable. Figures 16(a) and (b) present the pres- 
sure mode shapes for both the 23.8- and 2.3-hertz modes. One open-loop and two 
closed-loop, full-state feedback cases are shown. A propulsion system schematic is 
included along the abscissa to aid in locating the spatial coordinates of the various 
pressures throughout the system. Five maximums occur for the 23.8-hertz mode 
shape (fig. 16(a)); all have about the same magnitude in the open-loop case. For full- 
state feedback, as control weighting is decreased, the pressure maximums along the 
feedlines decrease relative to the highest maximum, which now occurs at the HPOP 
inlet location. Because each case is normalized with respect to its highest maximum, 
\no comparisons can be made between cases. But we can conclude that the control 
changes the way in which pressure perturbations affect the pressure distribution within 
the system. An increase in the amount of control decreases the effect of a perturba- 
tion on HPOP inlet pressure, as discussed in the section Frequency- Domain Behavior 
of Modal Control Designs. 

The 2.3-hertz mode shape (fig. 16(b)) was less sensitive to control than was the 
23. 8-hertz mode. However, the modal damping ratio for an r 2 of 0. 1 was almost 
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1000 times greater than the open- loop ratio . 

The flow mode shapes in figures 16(c) and (d) are like those in figures 16(a) and 
(b). Generally, flow is 90° out of phase with pressure, as would be expected. Again, 
increased control tends to alter the mode shape away from the open-loop, organ-pipe- 
like pattern. No clear pattern exists for the modal amplitude of control flow (state 29) 
as control weighting is decreased. This is true for both the 23.8- and 2 . 3-hertz 
modes . 


Sensitivity of Full-State Feedback Designs 

The modal design approach has been shown to be a very selective method. In it, 
the control is required to modify the system frequency response only in regions about 
the (two) selected modes. Because of this selectivity, the resulting designs could be 
overly sensitive to parametric variations and large increases in modal damping might 
not be achieved if system parameters changed. Thus, eigenvalue calculations were 
performed for off-nominal cases, where the most critical system parameters were 
varied. 

For the pogo model of this study, there are a number of parameters whose values 
are imprecisely known. But the most critical ones, it is generally believed, are the 
structural modal frequencies. Thus, off-nominal runs were made in which the 2.3- 
and 23.8-hertz modal frequencies were varied to see how they influenced the closed- 
loop modal damping ratios. Gains for the lower gain modal design (for r 2 = 1. 0) were 
used throughout. Figures 17(a) and (b) show how the damping ratios for all modes less 
than 30 hertz varied as structural modal frequencies were varied ±10 percent from 
nominal. Modes greater than 30 hertz were influenced very minimally. In figure 17(a), 
the 2. 3-hertz mode was varied. The only closed-loop mode whose damping ratio varied 
appreciably was the (nominally) 2.89-hertz mode. However, its damping is well above 
16 percent over the range of variation. 

In figure 17(b) the 23.8-hertz modal frequency was varied ±10 percent. As a re- 
sult a number of damping ratios changed, the most troublesome being the 26. 3-hertz 
mode. However, its lowest damping ratio is about 1 percent, which is still an accept- 
able value. The 23. 8-hertz modal frequency appears to be a much more influential 
parameter than the 2.3-hertz mode. Thus, for any final design, it would be advisable 
to more thoroughly investigate the sensitivity of any designs to the 23. 8-hertz modal 
frequency . 
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Partial-State Feedback 


Although the results and modal designs just discussed are useful as baseline in- 
formation, such designs would be impractical to implement in full- state feedback form. 
Implementing a Kalman filter would be one solution to the problem . But, because of 
the rather selective nature of the performance index, we believe an adequate job could 
be done more simply by using partial-state feedback. To see what gains might be 
eliminated, control gains for designs with of 1 and 0. 1 were tabulated (table V). 

In scanning the table, it is obvious that the largest gains occur for states 30, 31, 40, 
and 41 . This is not surprising, as these states are associated with the 2. 3- and 23 . 8- 
hertz modes; that is, they are the modal displacements and velocities for both modes. 
These two modes were the only ones weighted in the modal performance index. Thus, 
the partial-state feedback designs considered for further analysis used only states 30, 
31, 40, and 41. 


Generation of Unmeasurable Modal State Variables 

Partial-state feedback involving only states 30, 31, 40, and 41 still cannot be im- 
plemented because these states are generalized structural displacements and velocities 
and are thus not directly measurable . The way to get around this problem is to look at 
the equation relating structural outputs (which are physical quantities and can be mea- 
sured) to the structural states. This equation is given in appendix B and is repeated 
here for convenience. 


y = C x + C x (32) 

J s s s s s K ' 

Matrix C„ accounts for the liquid-oxygen-tank outflow effects. Its only nonzero 
s 

elements are ^-(tb), i = 1, . . .,6, for the six-structural-mode model. However, 

1 _ o 

for the model considered, all ^(tb) are zero except for £^(tb), which is -2.16x10 
System frequency responses and eigenvalues were calculated by assuming that £^(tb) 
is zero, that is, by neglecting all tank outflow effects. Very little difference was noted 
with and without ^(tb); therefore it was assumed that 

y s = C s x (33) 

Now, from the C matrix (table VI), it is evident that not all structural outputs are 

b 

independent. In fact, only six independent structural outputs can be defined. This 
six-element output vector, denoted as y" , can be defined as 
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( 34 ) 


By using C as defined in appendix B and the modal data given in table II, the output 

1^1 ® 

vector y can be related to the six generalized structural velocities as follows: 
s 



- 1.267 

- 0.198 

0.054 

0.022 

0.022 

0.022 ' 

- 0.8 

- 0.245 

- 0.0787 

- 0.752 

0.370 

0.0134 

- 2.38 

0.239 

- 0.556 

- 1.71 

- 2.13 

2.63 

- 0.137 

- 0.334 

- 0.259 

- 7.74 

4.44 

5.34 

- 0.316 

- 0.352 

0.431 

- 0.372 

0.017 

- 0.09 

- 2.68 

0.50 

- 1.95 

- 0.12 

0.99 

- 1.56 
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( 35 ) 


or, more compactly, 


^s - C a x a 


( 36 ) 


where 


a 


( q r q2.q3.Q4.Q5> q 6 ) 


( 36 ) 


Next, so that velocities x can be reconstructed, C must have an inverse, which it 

a ’ a 

does. In fact, it was calculated to be 


- 0.7248 

1.852 

- 0.06638 

- 0.09701 

- 1.347 

- 0.3606 

0.09167 

- 13.70 

0.3960 

0.7353 

9.771 

2.505 

1.6650 

- 7.951 

0.1027 

0.4207 

6.506 

1.179 

1.119 

1.707 

- 0 . 2015 

- 0.1130 

- 2.621 

- 0.5447 

0.8755 

- 0.6522 

- 0.2674 

0.1380 

- 0.2871 

0 . 045 02 

0.9128 

1.821 

- 0.04165 

- 0.02731 

- 2.668 

- 0.6223 


( 37 ) 


23 



Thus, 


x a = V. < 38 > 

Finally, the two generalized velocities (x sl = q 1 and x 4;J = q^), which need to be re- 
constructed, can be obtained as 

^1 


% 

Therefore to generate or reconstruct each generalized velocity requires the sens- 
ing of six structural velocities and the weighted summing of these velocities where 

A 

the appropriate row elements of matrix C Q are the weights . The two generalized dis- 
placements can be generated by integrating the respective generalized velocities. In 
an actual application, the only available structural measurements would be from ac- 
celerometers. These accelerometer outputs would need to be high-pass filtered to re- 
move any low-frequency component, linearly combined according to equation (39), and 
then bandpass filtered at the associated modal frequencies to determine the recon- 
structed modal velocities. This process would eliminate concern over the effect of 
any modes present but not considered in the pogo model, except for those that might 
occur'in a filter passband. If an extraneous mode does occur in a passband, additional 
information would be required to separate out its effects. 

Other methods of generating the generalized velocities from the direct use of par- 
allel bandpass-filtered data from one or more structural accelerometers can be con- 
ceived of. The important aspect is that the information required about the generalized 
modes that can destabilize through pogo coupling be supplied to the controller. Then 
with this information, the control energy can be concentrated at the offending modal 
frequency. 

At this point, the partial-feedback modal designs could have been further simpli- 
fied. However, because it was not obvious which elements of the matrix could be 
dropped, this aspect was not pursued. In the following sections, modal designs are 
further evaluated on the basis of eigenvalue-damping- sensitivity mode shapes and are 
finally compared with minimum-rms designs. 


KJ 4. 
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Mode-Shape Analysis of Partial-State Feedback Design- 

The partial- state feedback design that was evaluated used reconstructed estimates 
of states 30, 31, 40, and 41 (q^, q^, qg, and q^). The full- state feedback modal de- 
sign obtained for a control weighting r ^ of 1 was used as the basis from which the 
modal partial- state feedback design was obtained. 

In the partial- state design, damping may be decreased because fewer states are 
fed back than in the full-state design. For instance, the damping on the 23. 8-hertz 
mode decreased from 3. 67 percent to 2.42 percent in going to partial- state feedback. 
However, damping on the 2.3-hertz mode increased slightly, from 4.49 percent to 
4.87 percent. At any rate, it can be concluded here that for the four states chosen for 
feedback, the dampings on the two modes of primary interest remain well above the 
allowed minimum of 1 percent . 

Flow and pressure mode shape plots like those in figure 16 are presented in fig- 
ure 18 to clarify what happens when a full-state feedback design is simplified. Here, 
full- and partial-state feedback mode shapes are compared for the 23.8- and 2.3-hertz 
modes. First, the 2. 3-hertz mode shapes for both flow and pressure are almost iden- 
tical for full- and partial- state feedback. From this and the eigenvalue damping num- 
bers, it can be concluded that little performance was lost in eliminating most of the 
gains insofar as the 2 . 3-hertz mode is concerned. 

Differences in shapes do occur, however, for the 23. 8-hertz mode. In figures 
18(a) and (b), for instance, the partial-state feedback mode shape looks like one for a 
lower gain design; that is, it is more like an open- loop shape. This agrees with the 
fact that the damping ratio is less for the partial- state case (again, more like an open- 
loop ratio). The flow mode shapes for the 23. 8-hertz mode show the most pronounced 
differences. Mode shapes in the vertical and horizontal feedlines are similar, with 
partial-state feedback giving a mode shape more like open loop. The main differences 
occur in the interpump duct for interpump duct flow Q g and suppressor flow Q & . The 
interpump duct flows have a lower modal amplitude with partial-state feedback, and 
again this control looks more like an open-loop control in this region. This fact can be 
observed by comparing the open-loop mode shape values for Q g from figures 16(c) 
and (d) with the value shown in figure 18 . 

Sensitivity of Partial-State Feedback Designs 

One of the drawbacks to partial-state feedback is increased sensitivity to paramet- 
ric variations. This is obvious if the modal damping sensitivity in figure 19(a) is com- 
pared with that of the fuH- state feedback case (fig. 17(a)) upon which the partial- state 
design is based. At nominal conditions, some damping ratios are less for the partial- 
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state feedback case, notably for the 2 . 89- and 23 . 3-hertz modes. But these modes are 
quite heavily damped and the rest of the damping ratios are almost identical for full- 
and partial-state feedback. The main deficiency in the partial-state case is that there 
is an appreciable decrease in damping of the 2 . 3-hertz mode (off nominal) . For full- 
state feedback, this modal damping is insensitive . However, at its worst, the damping 
ratio is still an acceptable 1.5 percent. 

The sensitivity of partial-state feedback to changes in the 23. 8-hertz mode is 
shown in figure 19(b) . Comparing it with figure 17(b) shows that except for the 8.6-, 
23.3-, and 2.89-hertz modes, damping ratios and sensitivities are approximately the 
same with the full- and partial- state feedback. Fortunately, the damping ratios for 
these three modes are all above 2 percent in the ±10-percent range of parametric vari- 
ation. Thus, partial-state feedback affords quite adequate closed-loop stability. 


Comparison of Modal Control and Minimum-rms Designs 

The main thrust of the modal control design approach has been, naturally, to judge 
all designs according to the positions of their eigenvalues without regard to input and 
output behavior. This contrasts with the minimum-rms approach where input and out- 
put behavior is of main concern. Both approaches are valid; but it is instructive to 
examine modal control designs in terms of their ms behavior. This was done for the 
full- and partial-state feedback modal designs; The rms response values were calcu- 
lated by using the same colored-noise disturbance that was used to evaluate minimum- 
rms designs. The results are shown in figure 20, where normalized rms values of P 7 
are plotted as a function of normalized suppressor flow rate. For comparison, the 
mi nimum -rms, full-state feedback designs are shown by the solid line. The modal 
full- and partial-state design results are shown as the circles and squares, respec- 
tively . 

A number of points can be noted in figure 20. First, the modal designs used rela- 
tively small amounts of suppressor flow. This correlates with the fact that actuator 
poles for these designs are moved very little from the open-loop locations. Second, 
these modal designs can be characterized as low-gain, stabilizing controls. Third, 
the rms performances of full- and partial- state modal designs are not radically differ- 
ent, especially for the lower gain designs. However, none of the modal designs are 
close in rms performance to the minimum-rms, full-state feedback designs. This is 
to be expected since the design criteria for modal control designs differ from those for 
minimum-rms designs. 
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Frequency-Domain Behavior of Modal Control Designs 

For further insight into the modal control designs, a frequency-domain analysis 
was made. Thus, modal designs can be compared with previously discussed minimum- 
rms designs based on transfer function representations. Figure 21 shows open- and 
closed-loop frequency responses (magnitude only) for a modal design with full-state 
feedback. The main result is that disturbance is suppressed only in frequency 

bands near the two critical modal frequencies. Contrasting this to the minimum- rms 
response shown in figure 7(a), where disturbance suppression occurs throughout the 
range of frequencies shown, again emphasizes the selectivity that is achieved with the 
modal design. This frequency selectivity is also apparent in figure 22, where the 
closed-loop response of the control to disturbance is shown. Only two peaks occur in 
this plot - near 2.3 and 23.8 hertz, a result that is in sharp contrast to the many- 
peaked curve in figure 7(b) for a minimum-rms design. 

Going from full- to partial- state feedback causes little change in the frequency 
response of P 7 to disturbance . The response for partial- state feedback (fig. 23) 
is similar to the response of the original full-state feedback design (fig. 21), a fact 
that suggests the rms performances would also be similar. This, too, is the case, as 
shown in figure 20, where the leftmost circle and square are the points for the cases 
in question. Figure 20 also shows that normalized suppressor flow is less for the 
partial-state feedback simplification than for full-state feedback. (Compare figs. 22 
and 24 and the areas under the frequency responses for commanded suppressor flow 
rate.) The plot for the partial-state feedback case (fig. 24) encloses less area than 
does the plot for the full-state feedback case. In addition, a number of peaks occur in 
figure 24 that do not occur in figure 22. Thus, for partial-state feedback, the control 
attempts to attenuate disturbances at frequencies other than those near the two critical 
modes. This type of behavior was not expected based on the performance index used. 
However, such undesired control activity is not of significant magnitude to be a prob- 
lem. 


Use of Kalman Filter with Modal Control Designs 

As a final variation on designs that could be implemented based on the modal con- 
trol concept, a Kalman-filter-based controller was investigated. For this study, it 
was assumed that only HPOP suction could be measured, as was done for the minimum- 
rms design evaluations. Modal control gains for an rg of 1 were combined with Kal- 
man filter gains for a ^ w /^ v of 195 to form a feedback compensator such as that 
shown schematically in figure 10. The frequency response for this compensator is 
shown in figure 25. 
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One significant feature of this compensator’s magnitude response (fig. 25(a)) is 
that only two major peaks exist. However, responses for minim um-rms designs (fig. 
11(a)) have many peaks in addition to those near 2.3 and 23.8 hertz. The fact that the 
response had many peaks probably caused difficulties when low-order transfer function 
fits were attempted. Because of the relatively simpler shape of the modal-control 
compensator's frequency response, low-order fits should be more feasible. The com- 
pensator's phase angle (fig. 25(b)) shows much less overall phase lag than do any of the 
phase angles in figure 11(b) . This fact also makes it easier to obtain satisfactory low- 
order transfer function fits . 

A study to obtain such reduced-order models was not attempted here. These re- 
sults do suggest, though, that a viable control design solution might be a single P ? 
measurement fed through a doubly resonant low-order compensator. 


CONCLUSIONS 

The feasibility of applying linear quadratic regulator techniques and Kalman filter- 
ing to the design of an active-pogo- suppression control system was studied. Two lin- 
ear quadratic, optimal control approaches were taken: One minimized the root-mean- 
square response of high-pressure-oxidizer-pump suction pressure to a random distur- 
bance. The other increased the damping on two system oscillatory modes that were 
deemed most critical to pogo oscillation. Ideal, full-state feedback designs were com- 
puted and then used as baseline cases. Physical reasoning was used to find practical, 
partial-state feedback designs. A number of potentially implementable, simplified 
controls were found: 

1. Partial-state feedback, minimum-rms designs: These designs require sensing 
of two flow rates and two pressures. They increase stability of most system modes 
and reduce the overall system response to disturbances. Also, sensor requirements 
are not severe except those for a liquid-oxygen flow sensor, which may be a problem. 
These control designs were less able to increase damping on the critical system modes 
than were the modal control designs. 

2. Partial-state feedback, modal control designs: Each design requires that a 
sufficient number of structural accelerations be measured and appropriately filtered 

in order to reconstruct desired modal state information. Implementation is thus some- 
what more complex than for the designs described by Seidel in NASA TM X-3368. Con- 
trol energy is concentrated in fairly narrow frequency bands, and large increases in 
damping are achieved for the two critical system modes. Parameter sensitivity is not 
a severe problem. 

3. Reduced-order-feedback compensator based on modal control design using a 
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Kalman filter: This design uses only a single pressure measurement and a doubly 
resonant feedback compensator. The feasibility of obtaining a suitable low -order com- 
pensator looks promising. 

In general, the linear quadratic regulator approach appears to be well suited to the 
design of active-pogo-suppression controllers. Further, it provides means by which 
the issue of control energy can be effectively dealt with. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 18, 1978, 

505-05. 
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SYMBOLS 

pogo system matrix, 43x43 
colored-noise matrix, 2x2 
propulsion system matrix, 29x29 

partition of A matrix, 29x29 
structural system matrix, 12x12 

partition of A matrix, 12x12 

o 

cross-sectional area, cm 
acoustic velocity, m/sec 
system input matrix, 43x1 
propulsion system input matrix, 29x1 

partition of A matrix, 29x12 
partition of A matrix, 12x29 

Q 

compliance, cm /N 

low-pressure-oxidizer-pump compliance, cm^/N 
high-pressure-oxidizer-pump compliance, cm^/N 
structural output matrix, 12 x 12 

structural derivative output matrix, 12x12 
structural velocity matrix 
discharge coefficient, dimensionless 
disturbance matrix, 43x1 
colored-noise disturbance matrix, 2x1 
propulsion system disturbance matrix, 29x1 
coupling matrix, 29x12 
coupling matrix, 6x29 
structural mass matrix, 12x6 
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F 


force, N 

Fq tank-pressure coefficient matrix, 6x12 

F generalized force matrix, 6x12 

jr • generalized force, N 

f mean mass flow of liquid oxygen, g/sec 

f element of F matrix 

4 

g gravitational constant, N/g 

H output matrix, 12x43 

H cn noise-coloring output matrix, 1x2 

Hp partition of H matrix, 12x29 

H partition of H matrix, 12x12 

i component of vector 

modal performance index 
J-^ mean- square performance index 

j component of vector 

K c control gain matrix 

K Kalman filter gain matrix 

C 

2 5 

L inertance, N sec /cm 

£ line length, cm 

M generalized mass, g 

m^+1 low-pressure-oxidizer-pump dynamic gain, dimensionless 

nig+l high-pressure-oxidizer-pump dynamic gain, dimensionless 

P pressure, N/cm 2 

2 

Py high-pressure-oxidizer-pump suction pressure variation, N/cm 

P? m measured value of P^ 

2 3 

modal pressure coefficient, N-sec /cm 

2 

p pressure in feedline lump, N/cm 

O 

Q volumetric flow rate, cm /sec 

q generalized displacement, cm 
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— * 3 

q volumetric flow rate in feedline lump, cm /sec 

5 

R resistance, N - sec/cm 

R^ weighting matrix on states 

R 2 weighting matrix on controls 

r 2 element of R matrix (scalar) 

T modal transformation matrix 

U intermediate matrix, 6x6 

u control vector 

u commanded suppressor flow rate, cm /sec 

V intermediate matrix, 12x12 

v white-measurement noise vector 

W ^ intermediate matrix, 12x12 

W 2 intermediate matrix, 12x29 

w white-noise disturbance vector 

w cn colored-noise disturbance vector 

x overall system state vector, 43x1 

x cn colored-noise state vector, 2x1 

x propulsion system state vector, 29x1 

P 

x g structural system state vector, 12x1 

x ff vector of generalized structural velocities, 6x1 

y output vector, 12x1 

y structural system output vector, 12x1 

s 

z structural displacement along horizontal axis, cm 

y weighting on modal state variable 

£ damping ratio 

tj modal state vector 

A block-diagonal matrix 

O 

p liquid-oxygen density, g/cm 

j rms value of ( ) 


32 



t c time constant, sec 

$ mode shape vector, dimensionless 

cp modal displacement coefficient, dimensionless 

jp Y power spectral density of v 

power spectral density of w 
co structural pole natural frequency, rad/sec 

Subscripts: 
a actuator 

bl bubble for low-pressure oxidizer pump 

b2 bubble for high-pressure oxidizer pump 

D7 disturbance at station 7 

d discharge line 

I inner pump duct 

i index 

J injector 

j index 

k index 

LI first (vertical) feedline 

L2 second (horizontal) feedline 

£ lump 

n structural mode 

PI low-pressure oxidizer pump 

P2 high-pressure oxidizer pump 

tb tank base 

tc thrust chamber and discharge line 

Superscripts: 
x x-direction 

z z-direction 

~ normalized value 

derivative with respect to time 
('■) element of inverse matrix 
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APPENDIX B 


STATE -VARIABLE SPACE SHUTTLE POGO MODEL 

This appendix describes the dynamic model of the space shuttle oxidizer (liquid 
oxygen) system used for this pogo controller design study. Although the model is 
basically the one described in references 2 and 4, it has been modified to meet the 
needs of this study as well as other studies (e.g. , ref. 3) conducted at Lewis. The 
model can be divided into a structural portion and a propulsion system portion. The 
structural model, which is a six-mode approximation describing the liquid-oxygen- 
feedline motion, is based on one developed by the Aerospace Corp. and uses the struc- 
tural mode shape data presented in references 2 and 4. The propulsion system portion, 
which was also developed by Aerospace Corp. but modified at Lewis, is a linear, 
lumped-parameter model of the liquid-oxygen feedlines, pumps, and thrust chamber 
dynamics . The two portions of the model are joined through drag and momentum force 
coupling terms . 


System Geometry 

For dynamic modeling purposes a simplified liquid-oxygen feedline-pump-chamber 
geometry, developed by the Aerospace Corp. , was used (fig. 26). The system consists 
of the liquid-oxygen tank connected to a 31.1-meter- (102-ft-) long, 0.25-meter- 
(0.82-ft-) diameter vertical feedline; an elbow; and then an 8. 47 -meter- (27.8-ft-) 
long, 0.25-meter- (0.82-ft-) diameter horizontal feedline connected to the low- 
pressure oxidizer pump (LPOP). A short interpump duct connects the LPOP with the 
high-pressure oxidizer pump (HPOP), which pumps the liquid oxygen through a dis- 
charge duct into the thrust chamber. The suppressor is located at the inlet to the 
HPOP. 

The Aerospace Corp. structural model defines seven liquid-oxygen- system struc- 
tural elements that can move independently of one another: the tank, the two feedline 
segments, the two elbows, the LPOP, and the interpump duct - HPOP - discharge 
line - thrust chamber combination. As shown by the velocity coordinates in figure 26, 
the tank and the two feedlines can each move in only their axial directions (x^, x^, 
and Zg), but the elbows and pump assemblies are free to move both horizontally and 
vertically. Propulsion system dynamics are characterized by volumetric flow rates 
Q| and pressures P^ throughout the system. The subscripts i correspond to station 
locations: tb, tank base; 1, first elbow; 2, entrance to horizontal feedline; 3, second 
elbow; 4, LPOP inlet; 5, LPOP exit; 7, HPOP inlet; and tc, thrust chamber. These 
station definitions are the same as those used in references 2 and 4. 
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Propulsion System Dynamic Equations 

The propulsion system portion is modeled as a lumped-parameter system . The 
vertical feedline is divided into nine lumps and the horizontal feedline into three lumps . 
Figure 27 shows a typical lump, or control volume, used in obtaining the feedline equa- 
tions. Shown is the i n lump, which occurs in the horizontal (second) feedline. The 
vertical feedline equations, which consider conservation of mass and momentum for 
each lump, can be written as 


R» 


h 


q _£l i 

4l T 1 

h h h 


R £ - p i , p i-l R £ A 1 . _ 0 o 

q. H x n i 1 _ 2, 3, 

Lq 1 Lq Lq U K 


9 


q 9 ^Ql 


Pi = ^( q i' q i+l) 


i - 1, 2, . . . , 8 


P 8= P ] 


P 9 ~ P 2 “ TT ( Q 1 Q 2 + A l x l + A l z l) 


Similarly, the horizontal feedline equations can be written as 


**£ - p 10 P 2 ^ A 2 • 


r *10 -T+f- + 


10 


£ 2 


R, 


1 - 
1 


P 11 P 10 , R £ A 2 • 
+ + Zq 


q ll “ “ T q ll T T T T T % 2 

Lq L £ Lg L £ 


- 4 • Rf _ p 12 P 11 R £ A 2 • 

q 12 “ Q 4 “ ~ q 12 + , T z £ 2 
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PlO " « ( q io q n) 

4 

p li = ( q ll ■ ^12) 

4 


The cross-sectional areas of the two lines are assumed to be equal (A^ = A 2 = A); 
thus the compliance per lump is 


c t =S*L 

pa 2 


where g is the gravitational constant, p is the liquid-oxygen density, a is the speed 
of sound in liquid oxygen, (=p£ /Ag) is the inertance per lump, R^ is the resistance 
per lump, and P tb is the tank-base pressure. 

The equation for the pressure at the inlet to the LPOP can be written by consider- 
ing the cavitation bubble compliance and the motions of the elbow and LPOP (station 
locations 3 and 4 in fig. 26) as they affect the fluid system. The result is 


Pl2 ~ P 4 “ ^4 “ ^ 



where is the cavitation bubble compliance and A^ is the LPOP exit area, which 

equals the interpump duct area. The equations for the interpump duct were written for 
the duct as one lump with resistance Rj, inertance Lj, and cavitation compliance C ^ 2 
at the HPOP inlet. The resulting equations are 


Q 5 - ^ + + 1)P 4 - P ? - (Rj + RpilQs + R p i A 3 z 4 + ^ 3 * 7 ] 


Pr = 


5 = ( m i + X ) P 4 - R pl Q 5 “ L pl Q 5 + R pl A 3 Z 4 
P 7 = 7T“ [ Q 5 “ 4c + Q a + Q d7 + < A 4 " A 3^l\ 

4)1 J 


where Lp 4 is the LPOP inertance, m^ + 1 is the LPOP gain, R pl is the LPOP re- 
sistance, Q a is the suppressor flow rate (at the HPOP inlet), is the flow distur- 
bance at the HPOP inlet, and A^ is the HPOP suction area. All parameters associ- 
ated with the pumps and the interpump duct were obtained from references 2 and 4 . 
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The equations for the discharge line and thrust chamber pressure P^ c and flow 
Q^ c , which complete the propulsion system model, can be written as 

^tc = T 7 7" 77 [ P 7 " P tc _ (R d + R J + R p2 )Q tc + A 4^7] 

L d + L J + S 2 

and 



where the following parameters were also obtained from references 2 and 4: Ly, the 
discharge line inertance; Lj, the injector inertance; L p2* the HPOP inertance; Ry, 
the discharge line resistance; Rj, the injector resistance; Rp2» HPOP resistance; 
R tc , the thrust chamber resistance; and t^ c , the thrust chamber time constant. 

Structural Dynamic Equations 

In the model, seven structural elements characterize the liquid-oxygen system. 
The primary structural output quantities of interest are the horizontal z and vertical 
x velocities of each element. The structural elements and their associated indexes 
are tb, tank base; £1, vertical feedline; 1, first elbow; £2, horizontal feedline; 3, 
second elbow; 4, LPOP; and 7, the interpump duct - HPOP - discharge line - thrust 
chamber combination. 

The structural model was developed as in references 2 and 4 except for the treat- 
ment of feedline drag forces. The modal equations can be written as 

M n(% + 2 *AA + < 4 %) - 3, n - 2 6 

where, for the n 1 mode, M n is the generalized mass, q n is the generalized dis- 
placement, £ n is the damping ratio, u> n is the natural frequency, and is the 
generalized force . 

The model in references 2 and 4 has six modes, which are assumed to adequately 
describe the system at a particular flight condition. In this model the structural ve- 
locities x^ and z ^ can be related to the six generalized velocities by 
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6 

*i = Z i=U,l,3,4,7,tb 

n=l 

5 

Zi = I 4 Z) (i)% i =12,1, 3,4,7 
n=l 


where <p^ x )(i) and <p^\i) are components of the mode shape vectors. For example, 
the n^ mode shape vector 4> n is written as 


$ 


n 


<4 X) (U) 

<4 x) d> 

^ X) (3) 

^ X) (4) 

„f ( tb) 

4 Z) (* 2 > 
*?>< 1 ) 
4 z ) ( 3 > 

^ Z, (4) 

4 Z) (7) 


where the ordering of structural velocities is arbitrary. The remaining structural 
output (tank pressure) can be related to the generalized accelerations, as shown in 
references 2 and 4, as 



6 


i=l 


where coefficients ^(tb) are given. 
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The generalized force on the n n mode was shown in references 3 and 5 to be 


^ + A 1 Z \<4 X) ( tb > + P tb A l« , n X) ' tb > + Z F i X> ^ X) (i) + Z ’(J) 


i =£1,1, 3,4,7, tb 

j =£2, 1,3, 4, 7 - ■ 

n = 1,2, . . .,6 

The first term accounts for forces generated by liquid-oxygen-tank outflow and the sec- 
ond accounts for forces generated by tank pressure variation. The terms and 

Fj z) are the actual forces on the structure due to flow rate, pressure, and structural 
velocities at the seven defined locations. These forces are made up of (1) drag forces 
on the feedlines, the interpump duct, and the discharge line; (2) momentum forces at 
the two feedline corners, the LPOP, and the HPOP; and (3) forces on the thrust cham- 
ber. 

In a manner like that used in references 2 and 4, these forces can be expressed in 
terms of flows, pressures, and structural velocities previously defined: 

(1) Vertical and horizontal feedlines: The modeling here differs from that in ref- 
erences 2 and 4 since the feedline models are lumped instead of distributed. The drag 
forces are 


F £ ? “ _A l R je( X q i + A 1 x £ 1 


(2) First elbow: 


F £ 2 } A 2 R £ q i " A 2 Z £ 2^ 

F i X) = - P 2 A 1 ' * ( 2 ^ + X l) 


I< Z) =-P 2 A 2-fh^-Zi 
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(3) Second elbow: 


F^ x) = P 4 A 2 + f ( 2 5^ + x. 


f! z ^ = P.A 0 + f | 2 — + z„ + 2x„ 


3 ~ 4 2 


J 3 ' 3 


(4) LPOP: 


F 4 X) * “ P 4 A 2 - f 1 2 — + x. 


F f = - P 5 A 3- f ( 2 ^- Z 4 
A 3 


(5) Interpump duct - HPOP - discharge line - thrust chamber: 

F 7 X) = \c c ?tc 

F 7 ) = A 3 R l( Q 7 “ A 3 z 7) + A 4 R d( Q tc ' A 4 Z 7) + P 7 A 3 " P c A 4 + 2 * ^ + A 4 P tc 

Here, f is the mean mass flow of the liquid oxygen, A^ c is the thrust chamber area, 
and Cj is the discharge coefficient. 


Suppressor Dynamics 


The suppressor piston actuator is modeled as a first-order lag process with the 
transfer function 

Q a (S) = u a (S) 

S + 0)„ 


where u a (s) is the commanded flow input that produces the desired output suppressor 
flow rate Q„(s). 

d 
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Expressing Model in State -Variable Form 


The linear quadratic regulator approach to pogo suppressor design requires that a 
state- variable model of the standard form 

x = Ax + Bu + Dw 

y = Hx 

be obtained for the liquid-oxygen system . To do this, first a 12th-order structural 
system state vector x g is defined as made up of modal displacements and velocities. 

X I = (q r ■• ■>%>%) 

Then a 29th-order propulsion system state vector x^ is defined as 
*p = (‘Ip Pf* q 3 ,p 3’ q 4 ,P 4’ < V p 5 ,q 6 ,p 6 ,q 7 ,I V 

^8*^8* ^1* p 10’ q ll ,p ll’ ^4* ^4*^5* ^7’^tc*^a) 


Note here that suppressor flow rate Q a is included as part of the propulsion system 
state vector. 

For convenience in forming a state-variable model, the following vectors are de- 
fined: 


T A 


v g = (^i> • • •» &§) 6th-order vector of generalized forces 


y s _ (^ei» x l> x 3’ x 4 ,x 7 ,x tb ,z £2 ,z r z 3’ Z 4 ,Z 7’ ^tb) 


12th-order vector of structural- system outputs 


Also, the system disturbance w is selected as flow-rate disturbance . (Tank 
pressure is considered to be a structural output.) With these definitions, state- 
variable equations can be written for the structural system as 


k = A x_ + E v 
s s s s s 


y„ = C 0 x 0 + C 0 i 
J s s s s s 


v 0 = F„y. + D oi x + F„x„ 
s s J s sp p q s 
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by using the equations appearing on pages 35 to 40. Matrices A g , E g , C g , C g , F 
D a „, and F n are given in table VI. 

bp 4 

For the propulsion system, state equations can be written as 

VWV + WV 


where matrices A , B , D , and D „ are given in table VII. 

P P P ps 

Because of the x c terms on the right sides of the structural state equations, the 

b 

equations are not yet in "standard" state-variable form. However, by a straightfor- 
ward but lengthy substitution process, the state-variable equations 


x„ = A x + B x 
s s s s p 


x = Ax + B x + B u + D w 
p P P P s p p 


■ Vp + V s 


can be obtained, where the following matrices are defined: 

A s = A s + EgU(F fi Wj + F q A s ) 

A P = V D ps W 2 
B s = E s U ( F S W 2 + D S p) 

Bp = D pa W l 

H p * ^sV 1 + F , E S ) D sp 
» B = V{C S + C,[l + E s F q (I - E s UF q) ] A s } 


and the following intermediate matrices are defined: 


»= (I-FqEgf 1 


V = [l - C S E S (I + F q E s U)F s ] 


-1 
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w 




+ C S (! 


+ EF +EFEUF 
s q s q s 




C s E 


sd + F q E g U)D 


sp 


Colored-Noise Disturbance 

A second-order, colored-noise process, with natural frequency a; and damping 

vll 

t , was chosen to model the disturbance w. The state equations for this coloring can 
cn 

be written as 




or 


x = A x + D w 
cn cn cn cn 


w_„ = x _ 
cn cn cn 


Here, w is considered to be a white-noise disturbance. 


Overall State Equations 

Finally, by defining a 43rd-order overall state vector and a 12th-order output 
vector as 



the desired overall state equation becomes 
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x = Ax + Bu + Dw 


y = Hx 

where 



Parameter Values 

The nominal parameters that characterize the propulsion system are given in 
table I. Except for line inertance and compliance and actuator frequency, all were ob- 
tained from references 2 and 4. 

Bubble compliances and aod pump gains m^ + 1 and mg + 1 are es- 
pecially dependent on flight condition. The values shown in table I define the ejected 
parameter ranges. Structural parameters required are mainly the modal vector data 
plus modal frequencies, masses, and damping ratios. Table II presents these data, 
which are adapted from data in references 2 and 4, for the two flight conditions of 
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interest - end burn and after solid- rocket separation. These two conditions were 
chosen because pre lim inary analysis showed their pogo instabilities to be most pro- 
nounced. 


Nominal Worst-Case Model 

A nominal worst-case model was developed to reflect the dynamics of the two most 
critical flight conditions - end burn and after solid-rocket separation - and to achieve 
an open-loop model of the more unstable condition. The worst-case model is basically 
the reference 2 model at end burn except that the end-burn structural mode at 2. 81 
hertz (El) was replaced by a less stable mode (Al) from the after- separation condition 
(table II) . The frequency of the fifth end-bum mode was shifted from 26 . 8 hertz to 
23.8 hertz, making the system open-loop unstable. Also, the LPOP and HPOP gains 
were set to the maximum values shown in table I, and the pump inlet bubble compli- 
ances were set to the minimum values shown in table I. 
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APPENDIX C 


LISTING OF SUBROUTINE POGMAT 


Subroutine POGMAT is written in FORTRAN and performs the necessary matrix 
manipulations to calculate matrices A, B, D, and H given, as data, the parameters 
listed in table I and the structural information contained in tables VI and VII. 


0300100 SUBROUTINE POGM AT ( I DBU 3) 

0.000200 3 THIS SUBROUTINE CALCULATES THE A, B, C, HH, AND D MATRICES FOB THE '4 1 3 F ORDER 

0000300 0 POGO PROBLEM. A LS 0, DI MS NS I ONS N, MM, NC, ND, NO, AND NU ARE INPUT FROM POGMAT 

0003400. C TO THE MAIN PROGRAM. 

0000500 IMPLICIT RE AL*3 (A-H, 0-5 , 6-7.) 

0000530. . .. R£AL*8 T, TAUC, TAUCR, TAUP, TEMPB, TEMPO, TS, TSPTR , TST, TT 

0000700 REAL* 9 L, ICONST, LL, NN, K C, LLR, LI, LP1 , L D, LJ, LP2, LILPR, 

0300900 1 KE, LDJP2, LPA 

0000900 COMMON /ABETC/ A(53,53), B(53,5), 0(50,50), 0(50,15), HH(5,50), 

0301000 1 QC (50,50) , NN (50, 5) , PCINV(5,5), 00(50,50), RRINV(5,5) 

0001100 COMMON /COM1/ ADBLE(2530), EX1(50,50), EX2(53,50), EX3(50,50), 

0001200 1 EX4 (50) , EXT ( 1 00, 1 0 3 ) , KC(2,50), A MBKC ( 50 , 50 ) , ANR(50), ANI(50), - 

0301300 2 EIGR (50) , EIGI(53), X( 100, 100), EIGCLP.(50), EIGCLI (50) , TS2(133),- 

000 1 4 0 C 3 AR(100), AI(100), XR(1 33, 1 00), TT(100, 100), A A A ( 1 00, 1 00) , 

0001530 4 CP (100) , Cl (100) , S(53,50), SS(50,50), SSS(53,50), FREO(5CO), 

0001600 5 AMP (500) , PHASE(500) , AMP1 (500) , PHA1 (500) , AMP2 (500), PHA2 (503) r - 

0001700 6 AHPSTP (1030) , PHA 5 T R ( 1 3 33) , AINV(50,50), AM B INV ( 5 0 , 50) , 

0301930 7 ABKHIN(50,50) , A3 X CEH ( 5 0 , 5 0) , IBL(IOO), IA(130), 18(130), 

0001 900 8 LEX(IOO), MEX (1 03 ) , KE ( 5 3 , 5 ) , PP (53 , 50) , XTTM E ( 1000) , U (1000) , - 

0001920 9 EX5(50,2) ,EX6 (5,50) 

C002000 COMMON/COM2/ AP(33,30), BP(30,2), DPS(30,20), CP(30,30), 

0002100 1 AS (20,20) , ES(20,10), 05(23,15), CS(20,20), FS(10,20), PP(10,2), - 

0002200 2 DSP ( 10, 30) , PE(13(, DP(30,15), 7,B(10) 

0002400 COMMON/COM3/ CUSX3 (2,53) , CUU(2,50), BU(50,53|, AU(50,50), 

0302500 1 CTOT ( 100, 100) , DT 3 T ( 1 30 , 1 5 ) , XOLD(IOO), XH?B(100|, TEMPC (2,5 0) , - 

0002600 2 TEMPB (50,50) , ATOT ( 1 00, 1 00) 

0002700 COMMON/DIMS/N, NM, NC, ND, NO, ND, N2 

0002300 COMMON /STROCT/ NS, NP 

0002900 NAMELIST/PARIN/ RP1 , CB1 , CB2, BSTAR, AREA, AREA1, AREA2, A RE A3 , - 

0003000 1 AREA4, L, RH3, G, EL, Rt, RD, RC, RJ, RP2 , LI, LP 1 , LD, LJ, LP2, - 

0003100 2 PHI, PM2, R3, R4, TAUC, PBAR, ATCF, N, NM, SU, NS, NSD, NSF, N P3, - 

0303200 3 NC, ND, CS, FE, ZE, E3, HA, SPDSND , I ACT , AP7 , I AS 

0003220 ; 

0003240 C ***************** DATA A****************************************************,,** 

0003260 C 

0003300 IAS = 0 

0003400 WRITE (2,2025) 

0003500 READ (2,4568) IAS 

0003600 WRITE (2,4533) 

0003700 READ (2,4568) TACT 

0003800 4568 FORMAT (II) 

0003900 55 CONTINUE 

0004000 R P 1 = .021424D0 

0004100 CB 1 = . 25D0 

0004200 CB2 = . 1 0DD 

0004300 AREA = 75.7D0 

0004400 AREA1 = 75.7D0 

0034500 AREA2 = 75.7D0 

0004600 L = 136. 0D0 

0004700 RHO = .03813D0 

0004800 G = 386. 4D0 

0004900 BSTAR = 2.434D-5 

0005000 RL = 9.85D-5 

0005100 WA = 314.16D0 

0005200 RC = . 103D0 

0005300 R I = . 004532D0 

030540C pn = .02C1RRD0 

0005500 RJ = . 146672D0 

0035600 RP2 = . 1 08769D0 
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0005700 

LI = .5356D-3 

0005300 

LP 1 = .4120-4 

0005900 

LD = .618D-3 

0006000 

LJ = .2060-3 

0006100 

LP2 = .1030-3 

0006200 

PM 1 = . 306D0 

0006300 

PM2 = . 42D0 

0006400 

T AUC = . 45D-3 

0006500 

ESALL = 1.0D0 / 

0006500 

R 4 = . 02018800 

0006700 

R 3 = . 004532D0 

0005300 

IF (IAS .NE. 1) 

0006900 

PM1 = .8000 

0007000 

PM2 = . 54D0 

0007100 

CB 1 = I.ODO 

0007200 

CB2 = . 1 10DO 

0007300 4500 

CONTINOE 

0007400 

N = 41 

0007500 

NM = 1 

0007600 

NS = 12 

0007700 

NSO = 12 

0007800 

NS F = 6 

0007900 

NP = 29 

0008000 

NPO = 29 

0008100 

NC = 1 

0003200 

ND = 1 

0008300 

NMAX = 50 

0003400 

NMMAX = 5 


GO 


0008500 

0008600 

0008700 

0008800 

0008900 

0009000 

0009100 

0009200 

0009300 

0009400 

0009500 

0009600 


NPMAX = 30 
NPOMAX = 30 
NSMAX = 20 
NSOMAX = 20 
NSFMAX = 10 
NCMAX = 5 
NDMAX = 15 
NOMAX = 50 
N2HAX = 100 
N2 = N ♦ N 
NO=NSO«-NPO 
NTOT-N + NU 


0009700 


DO 

4502 

J 

= 1.N 

0009800 


DO 

4502 

I 

= 1,NM 

0009900 

4502 . 

HH 

(I,J) 

= 

0. ODO 

0010000 


DO 

2001 

J 

= 1,NS 

0010100 


DO 

2001 

I 

= 1,NSO 

0010200 

200 1 

CS 

(I. J) 

= 

0. ODO 

0010300 


DO 

2002 

J 

= 1,NSO 

0010400 


DO 

2002 

I 

= 1 , N P 

CO 10500 

2002 

DPS (I , J) 

= 

= 0 . 0 DO 

0010600 


DO 

2003 

J 

= 1 , NSF 

0010700 


DO 

2003 

I 

= 1 , NS 

0010800 

2003 

ES (I, J) 

= 

0. ODO 

0010900 


DO 

2004 

J 

= 1,NFO 

0011000 


DO 

2004 

I 

= 1.NSF 

0011100 

2004 

DSP (I, J) 

= 

= 0.0D0 

0011200 


DO 

2005 

J 

= 1 ,NP 

0011300 


DO 

2005 

I 

= 1,NPO 

0011400 

2005 

CP 

<I,J) 

= 

0. ODO 

0011500 


DO 

2006 

J 

= 1,NS 

0011600 


DO 

2006 

I 

= 1 ,NS 

0011700 

2006 

AS 

(I, J) 

= 

0. ODO 

0011800 


DO 

2007 

J 

= 1,NSO 

0011900 


DO 

2007 

I 

= 1 r N S F 

0012000 

2007 

FS 

(Ir J) 

= 

0. ODO 

0012100 


DO 

2008 

J 

= 1,NP 

0012200 


DO 

2008 

I 

= 1 ,NP 

0012300 

2008 

AP 

(Ir J) 

= 

0. ODO 

0012400 


DO 

2009 

J 

= 1 ,NC 

0012500 


DO 

2009 

I 

= 1,NSF 

0012600 

2009 

FP 

(I,J) 

= 

0. ODO 

0012700 


DO 

2021 

J 

= 1,NC 

0012800 


DO 

2021 

I 

= 1 ,NP 

0012900 

2021 


= 

0. ODO 

0013000 


DO 

2022 

J 

= 1 , ND 


. 3 DO 
TO 4500 



0013100 


DO 2022 

I = 1,NP 

0013200 

2022 

DP(I,J) 

= 0. 0 DO 

0013300 


DO 2023 

J = 1 , H D 

0013400 


DO 2023 

I = 1 ,NS 

0013500 

2023 

DS (X, J) 

= 0.0D0 

0013600 


CS (1 , 2) 

= - . 156D3 

0013700 


CS (2, 2) 

= - . 33D0 

0013800 


CS (3,2) 

= - . 33D0 

0013900 


CS ( 4, 2) 

= . 255D0 

0014000 


CS (5, 2) 

= . 255D0 

0014100 


CS (6,2) 

= - .209DO 

0014200 


CS (7 ,2) 

= .0367D0 

0014300 


CS ( 8, 2) 

= .0367D0 

0014400 


CS (9 , 2) 

= .0367D0 

0014500 


CS ( 10 , 2) 

= . 846D0 

0014600 


CS (11 ,2) 

= . 846D0 

0014700 


IF (IAS 

• NE. 1) G3 TO 

0014800 


CS(1,2) 

= . 529D0 

0014900 


CS (2, 2) 

= . 208D0 

0015000 


CS (3,2) 

= . 209 DO 

0015100 


CS ( 4, 2) 

= 1 . 23 DO 

0015200 


CS (5,2) 

= 1.28D0 

0015300 


CS (6, 2) 

= . 257D0 

0015400 


CS (7,2) 

= .09400 

0015500 


CS (8, 2) 

= . 094D0 

0015600 


CS (9, 2) 

= . 094 DO 

0015700 


CS (10,2) 

= 1.68D0 

0015800 


CS ( 1 1 ,2) 

= 1.68D0 

0015900 


CS (12, 2) 

= - 2. 160-3 

0016000 

2026 

CONTINOE 

0016100 


CS ( 1, 4) 

= - . 198D0 

0016200 


CS (2,4) 

= - .24503 

0016300 


CS ( 3, 4) 

= - . 245DD 

0016400 


CS (4,4) 

= . 2 39 DO 

0016500 


CS (5, 4) 

= . 239D0 

0016600 


CS (6,4) 

= - . 334D3 

0016700 


CS (7, 4) 

= - . 352D0 

0016800 


CS (8,4) 

= - . 352D3 

0015900 


CS (9, 4) 

= - . 352D0 

0017000 


CS ( 1 0 , 4 ) 

= . 50D0 

0017100 


CS ( 1 1 , 4) 

= . 50D0 

0017200 


CS ( 1 , 6) 

= • 054 DO 

0017300 


CS (2,6) 

= - .078700 

0017400 


CS (3,6) 

= - .078703 

0017500 


CS (4, 6) 

= - . 556D0 

0017600 


CS (5,6) 

= - . 556D0 

0017700 


CS (6.6) 

= - . 259D0 

0017800 


CS (7,6) 

= . 431D0 

0017900 


CS (8, 6) 

= . 43 1 DO 

0018000 


C S (9 , 6 ) 

= .43100 

0013100 


CS (10,6) 

= - 1 . 9500 

0018200 


CS (11 ,6) 

= - 1.9503 

0019300 


CS ( 1, 8) 

= . 022D0 

0018400 


CS(2,8) 

= - . 752D3 

0019500 


CS (3,8) 

= - . 752D0 

0018600 


CS (0,8) 

= - 1.7 IDO 

0019700 


CS (5, 8) 

= - 1 . 7 IDO 

0018800 


CS (6,8) 

= - 7. 7 4 DO 

0019900 


CS (7, 8) 

= - . 37200 

0019000 


CS (8 , 3) 

= - . 37203 

0013100 


CS (9, 8) 

= - . 37 2D0 

0019200 


CS (10 ,8) 

= - . 1 2 DO 

0013300 


CS (11,8) 

= - . 12D0 

0019400 


CS (1 , 10) 

= .22200 

0013500 


CS (2, 10) 

= .3700 

0019600 


CS (3,10) 

= . 37D0 

0019700 


CS (4, 10) 

= - 2. 1300 

0019800 


CS (5,10) 

= - 2.13D0 

0019900 


CS (6, 10) 

= 4.443D0 

0020000 


CS (7,10) 

= . 017 DO 

0020100 


CS (8, 10) 

= . 0 1 7D0 

0020200 


CS (9,10) 

= . 0 1 7D0 

0020300 


CS (10, 10) 

= . 99D0 

0020400 


CS (11 , 10] 

i = . 9 9 DO 
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0020500 


CS(1, 12) 

= .02200 

0020600 


CS ( 2 , 1 2) 

= .013403 

0020700 


cs ( 3, 12) 

= .013400 

0020800 


CS (4,12) 

= 2.63D0 

0020900 


CS (5, 12) 

= 2.63D0 

0021000 


CS (6,12) 

= 5.336D0 

0021100 


CS(7, 12) 

= - .0900 

0021200 


CS (8 , 1 2) 

= - . 09 D3 

0021300 


CS (9, 12) 

= - . 09D0 

0021400 


CS (10,12) = - 1.56D0 

0021500 


CS ( 11 , 12) = - 1 . 5500 

0021600 


F8(1) = 

2.8100 

0021700 


FE(2) = 

4 . 7 3 DO 

0021900 


FE (3) = 

8.4900 

0021900 


FE (4) = 

22. 5 DO 

0022000 


FE(5) = 

26.800 

0022100 


FE(6) = 

24.003 

0322200 


ZE(1) = 

.0100 

0022300 


7,E(2) = 

.0100 

0322400 


ZE(3) = 

.0100 

0022500 


ZE (4 ) = 

.0100 

0022600 


ZE(5) = 

• 01D0 

0022700 


ZE (6) = 

.0100 

0022800 


FBAR = 2 

. 29 DO , . ’ 

0022900 


ATCF = 158.000 

0023000 


AREA3 = 

31.200 

0023100 


AREA4 = 

12.6D0 

0023200 


AREAP = 

AREA 3 ... 

0023300 


A P7 = - 

.1480-2 

0023400 

4589 

CONTINUE 


0023500 


IP (IAS 

.NE. 1) GO TO 5)00 

0023600 


FE ( 1) = 

2. 300 

0023700 

5000 

CONTINUE 


0323720 

c 



0023740 



0023760 C 






0023800 

CL = BSTAR * AREA * L 




0023900 

LL = RHO * L / G / AREA 




0024000 

ROL = 

RL / LL 




0024100 

CLR = 

1.0D0 / CL 




0024200 

LLR = 

1.0D0 / LL 




0024300 

RAL = 

RL * AREA / LL 




0024400 

CB1R = 

1 . 0D0 / CB1 




0024500 

CB2R = 

1.0D0 / CB2 




0024600 

LILPR 

= 1.0D0 / (LI ♦ LP1) 




0324700 

LDJP2 

= LD + LJ ♦ LP2 




0024800 

TAUCR 

= 1.0D0 / TAUC 




0324900 

A 6 = (RD + RJ + RP2) / LDJP2 




0025000 

ARM1 = 

AREAl » RL 




0025100 

ARN2 = 

AREA2 * RL 




0025200 

AR3 = 

ARE A3 * R3 




0325300 

A32R = 

AR3 * AREA3 




0025400 

A 1 2R = 

AREAl ♦ AREAl * RL 




0025500 

A22R = 

A REA 2 * AREA 2 * RL 




0025600 

FBA 1 = 

2.000 * FBAR / AREAl 




0025700 

FBA 2 = 

2.000 * PBAR / AREA2 




0025800 

FBA3 = 

2.0D0 * FBAR / AREA3 




0025900 

FBA4 = 

2.0D0 * FBAR / AREA4 




0026000 

PdLIB 

= (1.0D0 ♦ PM) * LI * 

LILPR * AREA3 


0326100 

RLNRL 

= ( B I * LP1 - RP1 * LI) 

' * 

LILPF. * AREA3 


0026200 

RLHRL 

= - RLMRL 




0325300 

LPA = 

LP 1 * LILPR * AREA3 




0026400 

APH2 = 

AREA4 * (PM2 *■ 1 . 3D)) 




0026500 

RA3L = 

RP1 * AREA3 * A ?. E A 3 * 

LI 

* LILPR 


0026600 

RLA3 = 

RI * LP 1 * A RE A 3 * AFEA3 

* LILPR 


0026700 

A26B = 

APM2 * (LJ / LDJP2 - 1 

.0D0 ♦ LP2 / LD3P2) ♦ AREA 3 

0025800 

A27B = 

AREA4 * (-LP2 * A5 ♦ RP2 

♦ RJ - LJ * AS ♦ R4) 

- FBA 4 

0026900 

A28B = 

AREA4 * (-LP2 / LDJP2 

♦ 1 

.0D0 - LJ / LDJP2) 


0027000 

FI OB = 

FBAR - RA3L 




0027100 

FI IB = 

AREAU * AREA4 * (LJ * 

A6 

- RJ ♦ LP2 * A6 - R4 

- FP2) 

0327200 

DO 1 J 

= 1.N 




0027300 

DO 101 

1=1 ,NC 




0027400 101 

8(0,1) 

= O.ODO 




0027500 

DO 102 

I = 1 , ND 
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0027600 102 

D(J,I) = O.ODO 

0027700 

DO 1 I = 1 ,N 

0027800 

A (1,0) = O.ODO 

0027900 

C(I,J) = O.ODO 

0029000 

IF (I .EQ. J) C ( I, 3) = 1.0D0 

0028100 1 

CONTINUE 

0028200 

BP (29, 1) = BA 

0028300 

IF (I ACT .NE. 2) 30 TO 3 

0028400 

BP (25 , 1) - - Vi * L I * LILPR 

0028500 8 

CONTINUE 

0028600 

D (28, 1) = PC / TAO: 

0028700 

DO 4 1=1,21,2 

0028800 

AP (1+ 1 , I) = CLR 

0028900 

AP (1*2 ,1*1) = ILR 

0029000 

AP (1, 1*1) = - LLP 

0029100 

AP(I*1,I*2) = - CLP 

0029200 4 

AP ( I, I) = - BOL 

0029300 

AP (23 , 23) = - BOL 

0029400 

AP (23, 24) = - LLB 

0029500 

AP (24 ,23) = CB1B 

0029600 

AP ( 24 , 25) = - CB1R 

0029700 

AP(25,24) = LILPB ♦ (Pill ♦ 1.0D0) 

0029800 

AP (25, 25) = - (JI t BP1 ) * LILPR 

0029900 

A P (25 ,26) = - LILPR 

0030000 

AP (26 , 25) = CB2P 

0030100 

AP (26 ,27) = - CB2B 

0030200 

AP (27 , 26) = (PN2 ♦ 1 .000) / LDJP2 

0030300 

AP(27,27) = - (PD *■ RJ *■ RP2) / LDJ 

0030400 

AP ( 27 , 28) = - 1.000 / LD3P2 

0030500 

AP(28,27) = RC * r AUCR 

0030600 

AP ( 28 , 28) = - TAUOR 

0030700 

AP (29 ,29) = - HA 

0030900 0 


0030900 

IF (I ACT .NE. 1) 30 TO 5 

0031000 

AP (24 , 29) = CB1R 

0031100 5 

IF (I ACT .NE. 2) 30 TO 5 

0031200 

AP (25 , 29) = - LILPR * ( R I - HA « Li; 

0031300 6 

IF (IACT .NE. 2 .AND. IACT . NE. 3) I 

0031400 

AP (26 , 29 ) = CB2R 

0031500 7 

CONTINUE 

0031500 C 


0031700 

DO 2010 J = 1 ,NP 

0031900 

DO 2010 I = 1 , N P 

0031900 2010 

A (I, J) = AP(I.J) 

0032000 

DO 21 10 J = 1 ,NC 

0032100 

DO 2110 I = 1 , N P 

0032200 2110 

B(I,J) = BP (I, J) 

0032300 

AP ( 29 , 26) = - 1.00*8 

0032400 

DO 3001 I = 1,17,2 

0032500 3001 

DPS (I, 1) = - RAL 

0032600 

DPS (18,2) = AREA * CLR 

0032700 

DPS (24,3) = - AREA * C31 R 

0032800 

DPS (24,4) = - DPS(24,3) 

0032900 

DPS (19,7) = RAL 

0033000 

DPS (21,7) = RAL 

0033100 

DPS (2 3,7) = PAL 

0033200 

DPS (18,8) = DPS (13, 2) 

0033300 

DPS (24 ,9) = DPS (24,3) 

0033400 

DPS (24, 10) = AREA3 * C31P. 

0033500 

DPS (25,10) = RP1 * AREA3 * LILPR 

0033600 

DPS (25,11) = RI ♦ AREA3 » LILPP 

0033700 

DPS (26,11) = (AREA4 - A3EA3) * C32R 

0033800 

DPS (27 ,11) = AREA4 * A5 

0033900 

DPS (28,11) = - RC * ARSA4 * TAUCR 

0034000 

DPS (1,12) = LLR 

0034100 

DO 3003 J = 1 , NS 

0034200 

JJ = NP ♦ J 

0034300 

DO 3003 I = 1 ,NP 

0034400 

SUN = O.ODO 

0034500 

DO 3002 K = 1 , NS 0 

0034600 3002 

son = son ♦ dps(i,R) * :s(K,J) 

0034700 3003 

a (I , j j) = sun 

0034800 

DO 3004 I = 2, NS, 2 

0034900 

J = I / 2 
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0035000 3004 ES(I,J) = ESALL 

CO 35 1 00 E S ( 1 2 , 6) = 2 . 2 DO / 292. 3 DO 

0035200 IF (IAS .SB. 1) GO TO 3100 

0035300 ES (2, 1 ) = 1 . ODO/129 2. ODO 

0035400 3100 CONTINUE 

0035500 DO 3005 I = 1.NSP 

0035600 ' K = 2 * I 

0C35700 3005 DSP(I.I) = - ARM1 « CS(1,K) 

0035S00 DO 3006 J = 3,17.2 

0035900 DO 3006 I = 1,NSF 

0036000 3006 DSP(I.J) = DSP(I.I) 

0036100 DO 3007 I = 1 , NS F 

0036200 K = 2 * I 

0036300 DSP (1,1) = DSP (1,1 ) ♦ C3(12,K) 

0036400 DSP (I, 17) = DSP(I,17) - FBA 1 * CS(2,K) 

0036500 DSP (1 , 18) = - AREA 1 * C3(2,K) - AREA2 * CS(8,!C) 

0036600 FS (I, 1 ) = - A 1 2R * CS(1,K) 

0036700 PS (1, 2) = - FBAR * CS(2,K) 

0036800 FS (1,3) = PBAP * C3(3,K) * 2.0D0 * FBAR * CS(9,K) 

0036900 FS (1, 4) = - FBAR * CS(4,K) 

0037000 F S (1,6) = AREA1 * CS ( 1 2 , K) 

0037100 FS(I,7) = - A22P * CS(7,K) 

0037200 FS (I, 8) = FBAR * :S(8,K) 

0037300 FS (1,9) = FBAR * C5(9,K) 

0037400 FS (I, 1 0) = F10B * CS(10,K) 

0037500 FS (I, 11) - RLA3 * CS(10,K) ♦ F11B * CS(11,K) 

0037600 3007 FS(I,12) = AREA1 » CS (5 , K) 

0037700 DO 3009 J = 19,23,2 

0037800 DO 3009 I = 1,NSF 

0037900 K = 2 * I 

0038000 3009 DSP(I.J) = ARH2 * C3(7,K) 

0039100 DO 3010 I = 1,NSF 

0038200 K = 2 * I 

0038300 DSP (I, 19) = DSP(I,19) - FBA 2 * CS(8,K) 

0038400 EX4 (I) = CS(3,K) - C5('4,K) ♦ CS(9,K) 

0033500 DSP (1, 23) = DSP(I,23) ♦ FBA 2 * EX4(I) 

0033600 DSP (1 , 24 ) = AREA2 * FX4(I) - PML1B * CS(10,K) 

0038700 DSP (1,25) = CS(10,K) « (RLMRI - FBA3) ♦ CS(11,K) * ( A R 3 ♦ FBA 3) 

0033300 DSP (1, 26) = - CS ( 1 0 , K) * LPA ♦ CS(11,K) * A25B 

0038900 DSP (1,27) = CS(11,K) * A27B 

0039000 3010 DSP (1,28) = CS(11,K) * A23B CS(5,K) » A TCF 

0039 100 IF (I ACT .NE. 2) GO TO 3050 

0039200 DO 3040 I = 1 , NS F 

0039300 K = I * 2 

0039400 DO 3035 J = 1 , NC 

0039500 3035 FP(I,J) = - I.P1 * LI » »A * LTI.PR * AREA3 * CS(1C,K) 

0C39600 3040 DSP(I,29) = (FBA3 *■ A R3 1 * CS(11,FC) - L P 1 *■ LILPR * (RI - WA * Lr|- 

0039700 1 * AREA3 * CS(10,K) 

0039800 3050 CONTINUE 

0039900 DO 3011 J = 1 ,NP 

0040000 DO 3011 I = 1,NP0 

0040100 CP (1,0) = O.ODO 

0040200 IF (I .EQ. J) CP(I,J) = I.ODO 

0040300 3011 CONTINUE 

0040400 DO 3013 J = 1,NP 

0040500 DO 3013 I = 1,NS 

0040600 II = NP ♦ I 

0040700 SUM = O.ODO 

0040800 DO 3012 K = 1,NP0 

0040900 DO 3012 IL = 1,NSP 

0041000 3012 sun = sun + ES(I,IL) * DSP(IL,K) * CP (K , J) 

0041 100 3013 A (II, J) = sun 
0041200 DO 3014 I = 2, NS, 2 

0041300 K = I / 2 

0041400 AS (I, I- 1 ) = - 39.47600 * FE ( K) * FE (K) 

0041500 AS (I- 1 ,T) = I.ODO 

0041600 3014 AS (I, I) = - 12. 566D0 • ZE(K) * FE(K) 

0041700 DO 3016 J = 1,15 

0041800 JJ = NP ♦ J 

0041900 DO 3016 I = 1,NS 

0042000 II = NP ♦ I 

0042100 sun = O.ODO 

0042200 DO 3015 K = 1,NS0 

0042300 DO 3015 IL = 1,NSF 
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0042400 3015 SOM = SUM ♦ ES(I,IL) » PS(IL,K) * DS(K,J) 

0042500 3016 A(II,JJ) = AS(I,J| ♦ 5JM 
0042500 C 

0042700 C ADD EFFECT OP PRESSURE ON ACTUATOR POSITION FOR CONTROL AT. 3 

0042800 C 

0042900 IF ( I ACT .ME. 3) 30 TO 3025 

0043000 DO 3022 1=25, N, 2 

0043100 3022 A(29,I)= A(29,I) ♦ AP7 * 13 SAP * A(26,I) 

0043200 3025 CONTINUE 

0043300 DO 3018 J = 1,NC 

0043400 DO 3018 I = 1 , NS 

0043500 II = NP * I 

0043600 SOM = O.ODO 

0043700 DO 3017 K = 1,NSF 

0043800 3017 SUM = SUM ♦ ES(I,K) * PP(K,J) 

0043900 3018 B(II,J) = SDH 
0044000 HH (1,26) = 1.0D0 

0044100 WRITE (6,PARIN) 

0044200 2025 FORMAT (IK, 'FOR APTEP SEPARATION MODE, T TP E 1') 

0044300 3019 FORMAT (IX, 'A = ') 

0044400 4580 FORMAT (IX, 'ENTER ACTUATOR LOCATION ..10R20R3') 

0044500 501 FORMAT (IX, *DSP = *) 

0044600 502 FORMAT (IX, • FS = • ) 

0044700 503 FORMAT (IX, *ES = •) 

0044800 504 FORMAT (IX, 'DPS = •) 

0044900 505 FORMAT (IX, 'B = *| 

0045000 506 FORMAT (IX, »D =-•) 

0045100 507 FORMAT (IX, 'AS = •) 

0045200 508 FORMAT (IX, *PP = ') 

0045300 509 FORMAT (IX, 'CS = •) 

0045400 510 FORMAT (IX, 'CP = •) 

0045500 511 FORMAT ( 1 X , ' C = ') 

0045600 512 FORMAT ( 1 X ,' HH = •) 

0045700 C 
0045800 C 

0045900 2 PRINT OUT THE A MATRIX 

0046000 C 
0046100 : 

0046200 618 FORMAT ( I 2} 

0046300 WRITE (6,3019) 

0046400 CALL MATPRT (A, N, N, N N AX) 

0046500 WRITE (6,505) 

0046600 CALL MATPRT (B, N, NC, N M AX) 

0046700 WRITE (6,506) 

0046800 CALL MATPRT (D, N, ND, NMAX) 

.0046900 WRITE (6,511) 

0047000’ . CALL M ATPRT ( C, NO, N , NOM 1 X ) 

0047100 WRITE (6 ,512) 

0047200 CALL M ATPRT ( HH, NM, N , NM 1 1 X) 

0047300 IF (IAMAT .EQ. 0) 30 TD 379 

0047400 WRITE (3) <(A(I,J). J = 1,1), I = 1,N) 

0047500 379 CONTINUE 

0047520 C 

0047540 ; IF IDBUG .GT. 0 ,PD30 SUBMATPICES ARE PRINTED OUT 

0047560 C 

0047600 IF (IDBUG) 39,39,38 


0047700 

38 

WRITE (6,501) 




0047800 

0047900 


CALL MATPRT (DSP 
WRITE (6,502) 

, NSF 

, NPO, NSFMAX 

0043000 
004 3 IOC 


CALL MATPRT (FS, 
WRITE (6,503) 

SSF, 

NSD 

, NSFMAX) 

0043200 

0048300 


CALL MATPRT (ES, 
WRITE (6,504) 

NS, 

NSF, 

NSMAX) 

0043400 

0043500 


CALL MATPRT (DPS 
WRITE (6,507) 

, NP, 

NSO 

, N PM A X) 

0048500 

0048700 


CALL MATPRT (AS, 
WRITE (6 ,508) 

NS, 

NS, 

NSMAX) 

0043800 

0048900 


CALL MATPRT (FP, 
WRITE (6,509) 

NSF, 

NC, 

NSFMAX) 

0049000 

0049100 


CALL MATPRT (CS, 
WRITE (6,510) 

NSO, 

NS, 

NSOMAX) 

0049200 

0049300 

0049400 

0049500 

39 

CALL MATPRT (CP, 

CONTINUE 

RETURN 

END 

N PD, 

NP, 

NPOMAX) 
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TABLE I. - PROPULSION SYSTEM PARAMETERS 


2 2 

First feedline cross-sectional area, A^, cm ; in 

487; 75.5 

2 2 

Second feedline cross-sectional area, Ag, cm ; in 

487; 75.5 

2 2 

LPOP exit area, A^, cm ; in 

201; 31.2 

2 2 

HPOP suction area, A^, cm ; in 

83.3; 12.6 

2 2 

Effective thrust chamber area, A^c^, cm (in ) 

1020; 158 

Acoustic velocity, a, m/sec; in/ sec 

5.183; 20 400 

Lump compliance, , cm®/N; in® /lbf 

5.94; 0.250 

HPOP compliance, C^, cm^/N; in® /lbf 

. . 5.94 - 23.8; 0.25 - 1 

5 5 

LPOP compliance, cm /N ; in /lbf 

. 2.38 - 2.98; 0.1 - 0.125 

Mean flow of liquid oxygen, f, N- sec/cm; lbf - sec/in 

4.01; 2.29 

Gravitational constant, g, cm/sec 2 ; in/see 2 

. 981.5; 386.4 

Line lump length, f , cm; in. 

345; 136 

Lump inertance, Lj , N- sec^/cm®; lbf-sec 2 /in® 

• • 7 . 46xl0 -6 ; 1.77 xlO ” 4 

2 5 2 5' 

Inner-pump duct inertance, Lp N-sec /cm ; lbf- sec /in . . 

. . 2.25x10"®; 5.36 xlO” 4 

LPOP inertance, L pl , Nsec 2 /cm®; lbf-sec 2 /in® 

. . 1.73x10”®; 4. 12x10”® 

HPOP inertance, L 2 . N-sec 2 /cm®; lbf-sec 2 /in® 

. 4. 33x10”®; 1. 03xl0 -4 

Discharge line inertance, L^, N sec /cm ; lbf- sec /in . . 

. . 2 . 60x10”®; 6. 18xl0" 4 

Injector inertance, Lj, N-sec /cm°; lbf- sec /in 

. . 8.67x10"®; 2.06xl0" 4 

LPOP gain, + 1 

1.31 - 2.20 

HPOP gain, m 2 + 1 

1.42 - 1.54 

5 c 

Lump resistance, R^, N- sec/cm ; lbf- sec/in 

. . 4.14x10”®; 9.85X10” 5 

Inner-pump duct resistance, Rj, N- sec/cm®; lbf- sec/in 5 . . 

. . 1. 91xl0 -4 ; 4. 53xl0 -3 

5 5 

HPOP resistance, R^gj N- sec/cm ; lbf- sec/in 

4. 58xl0 -3 ; 0.109 

Discharge line resistance, R^; N- sec/cm®; lbf- sec/in® . . . 

. . 8 . 49xl0” 4 ; 2. 029xl0 -2 

5 5 

Injector resistance, Rj, N- sec/cm ; lbf- sec/in 

..... 6. 17 xlO 3 ; 0.147 

Thrust chamber resistance, R c ; N- sec/cm 5 ; lbf - sec/in® . . 

4. 33xl0" 3 ; 0.103 

3 3 

Liquid-oxygen density, p, g/cm ; lb/in 

.... 1.055; 3. 813xl0” 2 

Thrust chamber time constant, t c » sec 

4. 5xl0 ” 4 

Actuator pole frequency, « , rad/sec 

a 

100 ir 




TABLE II. - STRUCTURAL MODE DATA FOR END-BURN AND 


AFTE R-SO LID-RO CKET-SE PARATION CONDITIONS 


Parameter 



End bum 



After 


l_ : 






separation 


Mode 


I 

2 

3 

4 

5 

6 

1 




Aerospace Corp. 

designation 



El 

E2 

E7 

E30 

E34 

E35 

A1 

Modal displacement: 








<p| x) («l) 

-0.156 

-0. 198 

0.054. 

0.022 

0.222 

0.022 

0.529 


-0.330 

-0.245 

-0.0787 

-0.752 

0.370 

0.0134 

0.208 

<pfm 

-0.330 

-0.245 

-0.0787 

-0.752 

0.370 

0.0134 

0.208 

<p{*W 

0.255 

0.239 

-0.556 

-1.71 

-2.13 

2.63 

1.28 


0.255 

0.239 

-0.556 

-1.71 

-2.13 

2.63 

1.28 

<p| X) (tb) 

-0.209 

-0.334 

-0.259 

-7.74 

4.443 

5.336 

0.257 

<pj Z) ( *2) 

0 . 0367 

-0.352 

0.431 

-0.372 

0.017 

-0.09 

0.094 


0.0367 

-0.352 

0.431 

-0.372 

0.017 

-0.09 

0.094 

^ Z) (3) 

0.0367 

-0,352 

0.431 

-0.372 

0.017 

-0.09 

0.094 

<pf\ 4) 

0.846 

0.50 

-1.95 

-0.12 

0.99 

-1.56 

1.68 

<?i Z) (7) 

0.846 

0.50 

-1.95 

-0.12 

0.99 

-1.56 

1.68 

Modal pressure coefficient, 

0 

0 

0 

0 

0 

0 

-2 . 16 xlO -3 

^(tb), N-sec 2 /cm 3 
Generalized mass, M^ 

292 

292 

292 

292 

292 

292 

1282 

Damping, 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

Structural pole natural 

2.81 

4.73 

8.99 

22.5 

26.8 

27.2 

2.30 

frequency, co^/2n, Hz 









TABLE III. - NUMBER OF RIGHT-HALF-PLANE ZEROS 


Control 

weighting, 

Vn 

Ratio of normalized disturbance to measurement 
noise, ih /ib 

2 

1.95 

19.5 

195 


Number of right-half-plane zeros 

2xl0" 3 

17 

19 

17 

2xl0' 4 

17 

19 

15 

2xl0 -5 

13 

— 

2 

4xl0" 6 

11 

6 

0 
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TABLE IV. - MODAL DAMPING RATIOS FOR MODAL CONTROL DESIGNS 


[Weighting on modal state variable, y j> 1 f° r both critical modes; full-state 
feedback.] 

(a) Open loop (b) Closed loop 


Control weighting, r 2 

1.0 

0.1 

Natural 

Damping 

Natural 

Damping 

frequency, 

ratio 

frequency. 

ratio 

Hz 


Hz 


*2.315 

0.4485xl0 _1 

*2.335 

0.1386 

2.889 

.1651 

2.889 

.1651 

4.716 

,9858xl0“ 2 

4.716 

,9858xl0“ 2 

8.381 

,4765xl0“ 2 

8.381 

. 4765 xlO' 2 

8.603 

,6648xl0 _1 

8.608 

. 6648 xlO -1 

14.08 

,4492xl0 _1 

14.08 

.4492xl0 _1 

*23.84 

. 3673xl0 -1 

*23.99 

.1152 

22.51 

. 8965 xlO' 2 

19.42 

.4127 XlO' 1 

19.42 

.4127 xlO -1 

22.51 

.8965 xlO' 2 

23.31 

.1223 

23.31 

.1223 

26.32 

.1886 xlO' 1 

26.32 

.1886 xlO' 1 

28.10 

. 1808 xlO -1 

28.10 

. 1863 xlO -1 

32.11 

. 1990 xlO' 1 

32.11 

. 1990 xlO' 1 

35.88 

.1274 xlO' 1 

35.89 

.1274 xlO' 1 

39.31 

.7499xl0' 2 

39.31 

. 7499xl0 -2 

42.26 

. 4410xl0' 2 

42.26 

. 4410 xlO -2 

44.63 

.2615 xlO" 2 

44.53 

,2615xl0' 2 

46.36 

.16 18 xlO' 2 

46.36 

. 1618 xlO" 2 

47.41 

. 1099xl0' 2 

47.41 

. 1099xl0" 2 

50.11 

1.000 

50.11 

1.000 

61.40 

1.000 

61.40 

1.000 

100.00 

.7071 

100.00 

.7071 

333.1 

1.000 

333.1 

1.000 


Natural 

frequency, 

Hz 

Damping 

ratio 

100.00 

0.7071 

*2.312 

.3836 xlO' 2 

2.889 

.1651 

4.716 

. 9858 xlO' 2 

8.381 

.4765 xlO -2 

8.603 

.6648X10' 1 

14.08 

.4492 xlO' 1 

19.42 

.4127 xlO' 1 

23.31 

.1223 

22.51 

. 8965 xlO -2 

*23.83 

.1493 xlO' 2 

26.32 

. 1886 xlO' 1 

28.10 

. 1808 xlO -1 

32.11 

. 1990 xlO' 1 

35.88 

. 1274X10" 1 

39.31 

. 7499xl0' 2 

42.26 

. 4410X10" 2 

44.63 

.2615X10' 2 

46.36 

.16 18 xlO' 2 

47.41 

. 1099xl0 -2 

50.11 

1.000 

61.40 

1.000 

333.1 

. 1.000 


a Critical mode . 





State State Control weighting, r 2 

number 1 




1.0 

0.1 

Control gain 

Flow rate 

1 

0.0318 

0.0871 


3 

0.00152 

0.00412 


5 

-0.0289 

-0.0793 


7 

-0.0299 

-0.0825 


9 

-0.000572 

-0.00224 


11 

0.0303 

0.0826 


13 

0.0323 

0.0886 


15 

0.00333 

0.00978 


17 

-0.0281 

-0.0767 


19 

-0.0306 

-0.0844 


21 

-0.00215 

-0.00681 


23 

0.0298 

-0.0813 


25 

0.0603 

0.174 


27 

0.0168 

0.0612 


29 

0.0406 

0. 125 

Pressure 

2 

2.14 . 

6.76 


4 

2.16 

6.85 


6 

0.0758 

0.248 


8 

-2.02 

-6.37 


10 

-2.00 

-6.33 


12 

0.134 

0.414 


14 

2.27 

7.18 


16 

2.30 

7.28 


18 

0.202 

0.685 


20 

-1. 87 

-5.91 


22 

-1.87 

-5.91 


24 

0.235 

0.729 


26 

1.29 

4.07 


28 

0.0448 

0.157 

Displacement 

30 

1 

1540 

5200 


32 

-33.6 

-108 


34 

89.3 

292 


36 

246 

801 


38 

1070 

3410 


40 

-14 600 

-44 200 

Velocity 

31 

216 

677 


33 

-1.27 

-3.75 


35 

4.09 

12.1 


37 

1.24 

3.57 


39 

3.04 

8.47 


41 

23.5 

89.0 





















TABLE VI. 

MATRICES A s , E s> C s> Cg, F s> AND D gp 

0 1 

1 -| 

I 

-( 2 ^nl) 2 - 45r£ l £ nl 

1 

1 


jo 1 1 


l 

j -(2^) 2 - 4 ^2 f n2 L 

0 

X "j 0 1 


| -( 27rf n 6 ) 2 - 4 ^6 f n6 


12x12 


0 



1/M, 

o — 1 

0 


1/^2 L 


0 


0 



1/M 6 


12x6 


*1 1 *2 
1 


d> 

6 

' 

0 | 0 



__T_ 

0 _ 


12x6 
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TABLE VI. - Continued. 


( F 4 xl 2 1-1,2 ^ 

( f s\ r- A fo< 4 x) < ei) 

( F s). =-M X) W 

( F s). =f<pj X) (3) + 2f>{ Z) (3) 

i, 3 

( F s). 4 =-M X) < 4 > 

Hs-w.," 

(A), . -‘fP’w 

1, 8 

( F s) 19 -M Z) < 3 > 


tot , »S : 


(Z) 


i, 10 


I l + L pl 


(4) 


( F s), 


-a 2 r -a 2 

a 3 K 3 A 4 


i. 11 

Wl.18 

KL 29 ‘- 1 ' 2 6: 


V ^'K ,VV *“! 


^ z ) (7) + Ei_ R.A 2 <pp)(4) 

*1 + V 


= 0 


( D sp) i>2 


-^(tb) 


fo ) = 0 j = 4, 6, 8, 10, 12, 14, 16, 20, and 22 

p i,j 

i, J Ej 

( D sp). • = j = 3. 5, 7, 9, 11, 13, and 15 

J 

( D sp\,l 7 = “ A A^ X)(£1) 'f *i X>(1) 
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TABLE VI. - Concluded. 


( D =p) lls * A 2 E l”l Z, “ 2 >-f ’’‘ Z>,1) 

5 2 

( D s P \ 21 = A ^^ Z)(l2) 

( D sp\ )23 = A 2 R 1 ^i Z) < t2 ) + ~ [^ X) ( 3 > - } < 4 > + ^ Z) ( 3 >] 


( D sp ) i>24 = A 2 [«f } < 3 > + 4 Z> ( 3 ) + 4 X) < 4 )] ' ^7 A 3 ^ Z> 

( D s P ) i>25 =~ [^ z) < 7 >- ^ z) < 4 >] + Wi z) < ? > - 


( D sp) ; 


^ 1,26 iL^Lp! 


^-]a 3 ^ z) ( 4 ) + 


[ - A 4 (m 2 + 1)1, + . + /-Jal-i. 

_ L d +L j + L p2 \ L d + L 3 + 


( t \ \ = a 

/®p2 L d " 2L i R d + R 1 L d 2\ n 2f 

^ s P'i,27 4 

V L d + VV ) 4 A 2_ 


<?| Z) ( 7 ) 


= A . 


^ SP ^ 2 8 _ “ 4 \L d +L j + Lp2 

( D sp ) ij29 = 0 


' d ' V) + A tc c f¥ >f)( 7 ) 



0 

1 

Hi | 1 j 

j H« 

fF n ) 

' q '6xl2 

1 

1 

| 

Hi ! ! ! 

H. 

1 

1 

1 

1 |if 

1 j 


1 

1 

0 

Hi i 4 : 

! (4.6. 


where 


( f q). . = A i [^(tb)^ X >(tb) + ^(tb)<pj x )(tb)] + ^i (tb) ^j (tb) 


























































Thrust chamber 



jOrbiter vehicle 


Solid-propellant 
rocket motor t2) 


Location of liquid-oxygen tank 


Station 
8 ~ 

7 ~ 
5 -v 

4 - 

3 - 



-High-pressure oxidizer 
pump (HPOP) 

Pogo suppressor location 

— Interpump duct 

- Low-pressure oxidizer 

n. . mn /I DO D \ 


Station 1 -? 


Liquid- 

~~ f 


oxygen 

tank 


rFirst feedline 



z 

,--Second feedline 
r Station 2 


Structural motion 


(a) Overall geometry. 


(b) Schematic of liquid-oxygen system. 


Figure 1. - Space shuttle system for pogo analysis. 


Control 



Figure 2. - State-variable model for pogo analysis. 
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Thrust chamber 



w 


0 



Plant noise, w Measurement noise, v 



Optima! control 
gain matrix 


Figure 6. - Linear quadratic regulator solution with Kalman filter in feedback loop. 
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Open loop 
Closed loop 
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Normalized rms H POP suction pressure. Op , percent 



Figure 9. - Normalized root-mean-square high-pressure-oxidizer-pump suction pressure 
as function of normalized rms suppressor flow - Kalman filter on measurement P 7 m . 
Minimum-rms design; rms value of disturbance Q gj, 10 percent. 


Disturbance flow 



Feedback compensator, G c (s) 


Figure 10. - System configuration showing Kalman filter and optimal 
control gains as feedback compensator. 
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Figure 11. -Fre 
tion for differs 


disturbance to 







Normalized generalized displacement gain Flow gain 







Figure 13. - Normalized root-mean-square high-pressure-oxidizer-pump 
suction pressure as function of normalized rms suppressor flow - 
partial-state feedback designs compared with full-state feedback designs. 
Minimum-rms design; no measurement noise-, rms value of disturb- 
ance Q d7 , 10 percent. 
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Normalized amplitude ratio, — - 



0 10 20 30 40 50 

Frequency. Hz 


Figure 14. - Normalized frequency-response amplitude ratio of output 
P 7 to disturbance (^7 - minimum-rms, partial-state feedback de- 
sign where control weighting r 2 is 2xl0"3. 



0 10 20 30 40 50 

Frequency, Hz 


Figure 15. - Normalized frequency-response amplitude ratio of output P 7 
to disturbance - where 23.8-hertz modal frequency has been in- 
creased by 10 percent. Minimum-rms, partial-state feedback design 
based on full-state feedback design where control weighting r ? is 
2x10-3. 
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Closed-loop modal damping ratio. ? 
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Figure 19. - Sensitivity of closed-loop modal damping ratios to critical structural mode frequency - square value of disturbance Q - 

partial-state feedback design based on full-state feedback design for control weighting r 2 of 1. 0. . 10 percent. D ’ 




Frequency. Hz 

Figure21. - Normalized frequency-response amplitude ratio of output P ? 
to disturbance Qqj - modal control. full -state feedback design. Control 
weighting, r ? 1.0. 



0 10 2D 30 40 50 

Frequency. Hz 


Figured - Normalized, closed-loop frequency-response amplitude 
ratio of control input u a to disturbance Q^j - modal control, 
full-state feedback design. Control weighting, r^ 1.0. 
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Normalized amplitude ratio, — - — P 7 (ioj) 

jO d7 (jw) Normalized amplitude ratio. - 



0 10 20 30 40 50 

Frequency, Hz 


Figure 23. - Normalized frequency-response amplitude ratio of output Pj 
to disturbance Qd7 - modal control, partial-state feedback design basea 
on full : 'state feedback design for control weighting r 2 of 1.0. 



0 1 0 20 30 40 50 

Frequency, Hz 


Figure 24. - Normalized, closed-loop frequency-response amplitude ratio 
of control input u a to disturbance Q d7 - modal control, partial-state 
feedback design based on full-state feedback design for control weighting 
r 2 of 1.0. 
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r- Oxidizer tank 



High-pressure oxidizer pump (HPOP)-^ 
Figure 26. - Schematic of liquid-oxygen-system model. 



Figure 27. - Feedline control volume. Lump length, [; 
pressure in i 1 * 1 lump (second feedline), p;; flow entering 
i th lump (second feedline), q;; horizontal feedline 
velocity, Z 12 ; cross-sectional area (second feedline), 
A 2 . 
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